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Abstract 

We obtain an explicit expression for the multipoint energy corre- 
lations of a non solvable two-dimensional Ising models with nearest 
neighbor ferromagnetic interactions plus a weak finite range interac- 
tion of strength A, in a scaling limit in which we send the lattice 
spacing to zero and the temperature to the critical one. Our anal- 
ysis is based on an exact mapping of the model into an interacting 
lattice fermionic theory, which generalizes the one originally used by 
Schultz, Mattis and Lieb for the nearest neighbor Ising model. The 
interacting model is then analyzed by a multiscale method first pro- 
posed by Pinson and Spencer. If the lattice spacing is finite, then 
the correlations cannot be computed in closed form: rather, they are 
expressed in terms of infinite, convergent, power series in A. In the 
scaling limit, these infinite expansions radically simplify and reduce 
to the limiting energy correlations of the integrable Ising model, up 
to a finite renormalization of the parameters. Explicit bounds on the 
speed of convergence to the scaling limit are derived. 
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1 Introduction 

1.1 The model and the main results 

The two-dimensional (2D) Ising model is an oversimplified description of a 
planar magnet in which the dipoles can point only in two directions and 
only the interactions among neighboring spins are considered. Remarkably, 
it can be solved explicitly, as first shown by Onsager [IB] : the free energy, the 
magnetization and several correlation functions in the absence of an external 
field can be computed in closed form [351 IMl EZl US] . 

From a physical point of view, there is no special reason to consider only 
nearest neighbor (n.n.) interactions: it is natural (and more realistic) to 
consider generalizations of the n.n. Ising model where the spins interact via 
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a finite range interaction of tlie form: 



Hm = -J ^ ^ O-xCTx+ae, ' Crx^^(x - y)cry = Hf^ + \Wm , (1-1) 

xeaZ|, i=l,2 {x,y} 

where J is a positive constant, C is a finite square box of side M witli 
periodic boundary conditions, a is tlie lattice spacing, cXx G {±1}, and are 
the two unit coordinate vectors on . The sum in the second term of Eq. fll.ip 
is over all unordered pairs of sites in := aZ|^; the interaction potential 
i;(x — y) is rotation invariant and has finite range, namely: f (x — y) = 0, 
V|x| > i?o := aMo, for a suitable positive integer Mq. The standard n.n. case 
corresponds to A = 0. 

Despite the fact that the models with Hamiltonian or Hm look "phys- 
ically equivalent" , Onsager's exact solution is crucially based on the assump- 
tion that A = 0: therefore, it is natural to ask how much of it survives the 
presence of the perturbation XWm- 

The thermodynamical properties of the system are obtained by averag- 
ing with respect to the Gibbs measure at inverse temperature /3: given an 
observable F{a) (here a = (cTxjxeAj), its statistical average is defined as 

where Qm = {±1}^^ is the spin configuration space. In particular, the av- 
erages (cTxi ■ ■ ■o':x.n)i3 L multipoint spin correlation functions. Taking 
L — )■ oo at {xi, . . . ,x„} fixed produces the infinite volume correlation func- 
tions: 

(f^xi ■■■cr^„)f^= lim (cTxi ■ ■ ■ fx^J^ ^ (1.3) 

the collection of which characterizes the infinite volume Gibbs state. Two key 
thermodynamic quantities are the free energy and the specific heat, defined 
as follows: the free energy per site in the thermodynamic limit is 

where Z = X]cref7j,f ^^p{~/^-^m(^)} the partition function. The specific 

heat is a = and can be obtained by summing over the energy 

density correlation functions: if we define the energy density operator as 
^xj := a""^o-x(Xx+ae,, then 

Cv = i(3JY Irni^^^ j rfxdy(£x,j;£yj')Ji , (l-^) 
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where dx is a shorthand for XlxeAj ^^"-^ ^'i indicates the truncated 

expectation: (exj;%j')Ji = {^^,j^y,j') p,L - ^^^,h) p,L(^yd2) p,L- Note that 
the thermodynamic functions hke the free energy and the specific heat are 
independent of the lattice spacing a while the correlation functions explicitly 
(if straightforwardly) depend on it. 

If A = 0, Onsager's exact solution shows that the free energy is analytic 
for 13 ^ /3c, where /3c = J"Hanh-^(y2 - 1) is the inverse critical tempera- 
ture. For /3 < /3c there is a unique infinite volume Gibbs state (the "high 
temperature state") with zero spontaneous magnetization (i.e. (cTx)^ = 0) 
and characterized by exponential decay of the multipoint spin correlations 
as the separation between spins is sent to infinity. For /3 > /3c there are ex- 
actly two pure Gibbs states (■)^, which generate all possible infinite volume 
Gibbs states by convex combinations (HED] (the "low temperature states"): 
they have a non zero spontaneous magnetization (i.e. (o"x)^ = ±m*(/3) ^ 0) 
and are characterized by exponential decay of the multipoint truncated spin 
correlations as the separation between spins is sent to infinity. At the crit- 
ical point, the rate of the exponential decay goes to zero proportionally to 
1/3 — /3c I and the correlation functions decay polynomially to zero, with spe- 
cific critical exponents] e.g., the spin-spin as (cr^ay) ~ (const.) |x — y| 
while the energy-energy correlation decays as {£x,j£y,j') ~ (const. )|x — y| 
asymptotically as |x — y| — )■ oo. Correspondingly, the specific heat diverges 
logarithmically ~ (const.) log |/3 — /3c| at the critical point. We remark 
that, even after having understood the integrability structure of the model, 
underlying the Onsager's computation of the free energy, the determination 
of the spin-spin critical exponent is very involved: ((TxCTq)^ is expressed as a 
determinant of a large matrix, of size increasing with the distance between 
spins, whose asymptotic behavior along special directions can be obtained 
via the use of Szego's lemma [13] or via the use of the analyticity structure 
of a set of exact quadratic difference equations that (o"xcro)^ satisfies exactly 
at the lattice level |44j . 

On the other hand, the computation of the multipoint energy correla- 
tions at A = is a relatively simple matter. Defining the truncated energy 
correlations as 

with 

2 

Z{A) = J2 exp{-/3HM{a) + / e^,,A^J, (1.7) 
it turns out that at A = the correlations (^xiji;--- ;^x™jv„)^2. 
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written in closed form for all finite a, L. Their expression further simplifies 
in the thermodynamic limit L — t- oo and in a scaling limit in which a and 
13 — (5c are simultaneously sent to zero: defining tc = ianh. jScJ = \/2 — 1, 
t = tanh/3(a) J and fixing /3(a) in such a way that 

t — tr aaia) 



with a{a) 7^ for a 7^ and lima^o <y{a) = m* G M, then, for all the m-tuples 
of non-coinciding points Xi, . . . , x^ G M^, m > 2, Xj 7^ x^ Vi 7^ j, 

where: (1) jk G {1,2} and the limit is independent of the choice of these 
labels; (2) n{i,...,n} is the set of permutations over {1, . . . ,n}, with n^n + 1) 
interpreted as equal to vr(l), and cok G {±}, with con+i interpreted as equal 
to ui, (3) S°(x) is given by 

^ J 27r k2 + (m*)2 V -^m* th - ' ^ ' 

which is a well known object in Quantum Field Theory (QFT): it verifies the 
Dirac equation in 1 + 1 dimensions, which describes a quantum relativistic 
fermion 

\^ -im* 01-102 J 

The derivation of Eq.f ll.8p will be reviewed below. It is based on a cor- 
respondence between the energy operator of the Ising model at the site x 
and the mass operator (z/7r)'?/'x,+V'x,- of a QFT of non-interacting Majo- 
rana fermions, with propagator ^{'ipyi,uj4'y,uj') = 0^w'(x "~ y)- The r.h.s. of 
Eq. fll.8p can be recognized as the multipoint truncated fermionic correla- 
tion (i/7r)™'E-^('?/'xi,+'^xi,-; ■ ■ ■ ; ■^x^.+V'x™,-), where E'^ is computed via the 
fermionic Wick rule and is given by the sum over the pairings ( "contractions" ) 
of the fields; each pairing is equal to the product of the propagators associ- 
ated with all the contracted pairs, times the sign of the permutation needed 
to bring the contracted fermionic fields next to each other. The outcome 
can be naturally represented in terms of Feynman diagrams and truncation 
means that only connected diagrams are considered; in graphical terms, the 
r.h.s. of Eq.f ll.8p can be thought of as a sum over simple loop graphs, see 

Figm 



5 




+ permutations 



Figure 1: Graphical representation of Eq. fll.22p in the case m = 6. The 
wavy hues represent the external fields ^xj, which can be thought of as 
being associated to the vertices. These are labeled by the space coordinate 
X of the pair of fields '?/'x,+'^x,- "exiting from the vertex". The fields are 
contracted two by two: the result of the contraction is represented by a 
solid line (labeled by the u, u' indices of the two contracted fields) and is 
associated with a propagator g°. 



Note that the thermodynamic limit L — )■ oo is performed at /3 = /3(a) 7^ 
Pc] after the thermodynamic limit /3(a) is sent to /3c and, depending on the 
speed at which /3(a) — )■ /3c, the resulting expression has a different effective 
mass m*, which is a finite real number. In particular, in the special case 
m* = 0, the propagator reduces to 

0L'(x)L_o = ^^. (1.11) 

'™ Xi + 10JX2 

In the literature, the explicit expression Eq.f ll.8p with m* = is usually re- 
ferred to as the critical multipoint energy correlation in the plane, see e.g. 
[T6] : the result is independent of the specific boundary conditions imposed 
on AJ. On the contrary, if the lattice spacing a is sent to directly at /3 = /3c, 
with L kept fixed, then one gets a different expression, which is sensitive to 
the specific boundary conditions imposed on A^: e.g. if periodic boundary 
conditions are considered, the multipoint energy correlations tend to the so- 
called critical correlations on the torus p!6] . 

All these results rely on the exact solution of the n.n. Ising model. As 
mentioned above, an exact expression for the energy correlations is known 
also at finite a and has a form similar to Eq.f ll.Sp . still expressed in terms of 
loop graphs but with a different (and more involved) lattice propagator; in 
the scaling limit, it simplifies and reduces to the nice expression Eq. fll.Sp . for- 
mally coinciding with the correlations of non interacting Majorana fermions. 
When A 7^ 0, the situation is radically different; there is no exact solution 
in that case, and in particular no explicit expression for the m-point energy 
correlations is known. If the temperature is well inside the high or low tem- 
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Figure 2: Schematic summary of the renormahzation of the correlation func- 
tions in the scahng hmit. The symbol ■ denotes a dressed vertex, propor- 
tional to a renormahzation factor Z{X). 

perature phase, one can pursue a perturbative approach based on cluster 
expansion methods, which allows one to prove quantitatively that the be- 
havior of the perturbed system is close to the exactly solvable one, see e.g. 
im HU [23]- The difficult and subtle case is when the system is close to or at 
the critical point. In this case, the cluster expansion argument breaks down: 
the perturbation theory is affected by infrared divergences related to the 
slow decay of correlations, which may a priori change the critical exponents 
and the nature of the critical point. Spencer |5l] and Pinson and Spencer 
|17] introduced a new point of view using a generalization of the mapping 
between 2D Ising models and fermionic systems, first discovered by Schultz, 
Mattis and Lieb |;52j. While the n.n. Ising model is equivalent to a system 
of non-interacting fermions, when A 7^ the equivalence is with a system of 
interacting fermions. By such mapping (which is reproduced and extended 
to more general interactions than those in |17] in Section H] below) one im- 
mediately gets a perturbative expansion for the energy correlations in terms 
of Feynman graphs of arbitrary order, see the first line of Figj2l contrary to 
the A = case, now graphs with any number of loops appear. 

The aim of this paper is to control and compute the scaling limit of these 
energy correlations: this means to control and compute the ground state 
correlations of the underlying interacting fermionic theory, uniformly in the 
effective mass m* . We show that the apparent infrared divergences at m* = 
can be resummed by rigorous Renormahzation Group methods: the outcome 



7 



is a renormalized expansion for the m-point energy correlations that is con- 
vergent for A small enough, as shown in [17] for the 2-point function. On 
top of it, and in addition to the estimates derived in fl7], we compute the 
explicit exact expression of the interacting multipoint energy correlations in 
the scaling limit and show that they coincide with an expression completely 
analogous to the one obtained for the integrable Ising model, up to a finite 
renormalization of the parameters. Graphically, the sum over infinite Feyn- 
man diagrams with an arbitrary number of loops in the first line of Figj2]can 
be shown to be equal to a sum over simple dressed loop diagrams, like the 
one in the second line of Figj2l up to subdominant corrections that vanish 
in the scaling limit. We stress that the use of fermionic methods is crucial 
to establish these results, which remained open problems for several years 
and appeared out of reach in terms of the original spin variables. Our main 
result is summarized in the following theorem. 

Theorem 1.1. There exists Aq > such that if \X\ < Aq the following is 
true. There exists an analytic function of X, called /3c{X), which is the critical 
temperature of the Ising model Eq. U.l]) . Defining tc(A) = tanh/3c(A)J, t = 
tanh/3(a)J and fixing /3(a) in such a way that 

t — tc(A) tc(A) a(T{a) 
tc(A) " tc(0) 2 ' 

with a{a) ^ for a ^ and lima_^o = m* G M, then, for all the m-tuples 
of non- coinciding points Xi, . . . , G M^, m > 2, Xj 7^ Xj Vi 7^ j, 

lim (^xi ji; ■ ■ ■ ; ^x™,i™) J(a),L = (1-12) 

where: (1) ji G {1,2} and the limit is independent of the choice of these 
labels; (2) rTji^ is the set of permutations over {!,... ,n}, with n^n + 1) 
interpreted as equal to tt{1), and Uk G {±}, with Un+i interpreted as equal to 
uji; (3) the dressed propagator in the scaling limit is 

g(x) = Z(A) / 7^ ,,, , ,U ' U_ ' (I'^S) 



271 k2 + [Z*(A)m*]2 \-tZ*{X)m* ih - k2 

where the renormalizations Z{X),Z*{X) G M are analytic functions of X, 
analytically close to 1. 
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Therefore, even in presence of an interaction, the continuum hmit of 
the energy operator can still be identified with the mass operator of a free 
fermion, up to a finite renormalization of the mass and of the wave func- 
tion. Note that the case m* = is included in the theorem. The critical 
temperature is renormalized by the interaction, too; on the contrary, the en- 
ergy critical exponents are protected against renormalization. The theorem 
is proved by combining the methods of Pinson and Spencer [171 with 
the analysis of the massive Thirring model developed by Benfatto, Falco and 
Mastropietro [8J. The proof of the existence of the scaling limit behind Theo- 
rem [LT] also allows us to fully control how fast the scaling limit of the energy 
correlations is reached, as stated explicitly in the next theorem. 

Theorem 1.2. Given A, /3(a) and Xi, . . . , as in the statement of Theorem 
defining 



s(<^)(x) = Z(A) 



dk e-'^"" / D+(k) iZ*{X)a{a 



__7V_ 7r-|2 
a ' a 



(1.14) 



we have: 



lim (exiji; ■ ■ ■ ; £x„,>„)L),l = (1-15) 



L— >oo 



2m \7i 



7r(fc) — X7r(fc + l)J 



+ 



1 JmJ ; 



where, denoting by the diameter of x. = {xi, . . . , x^} and by 6^ the min- 
imal distance among the points in x, for all 9 G (0, 1) and e G (0, 1/2) and a 
suitable Ce^ > 0, the correction term R^""^ can be bounded as 



1 f a \ / o-v 



|/2('^)(xi,ji;---;x^,j„)|<C,>!— f-j i^) . (1.16) 



Note that the contributions in the second line of Eq. fll.15p tend to Eq. fll.l2p 
in the limit as a — )■ 0, in a way controlled by the difference ||g'^°)(x) — g(x)||, 
which is bounded by |x|~-'^e~('^°°'^*'^™*l^l (a|x|~^ + \a{a) — m* \ ■ |x|) as a — )■ 
and a{a) — )■ m*. Note also that the combinatorial factor ml in the r.h.s. 
of Eq. fl 1.161) is the optimal one and corresponds to the same combinatorial 
growth as m — >■ oo as the dominant term in the second line of Eq. fll.lSp . 

Moreover, let us remark that the proof of Theorem 11.21 tells more than 
what is stated above. In particular, it allows us to distinguish, among the 
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contributions to R^"'\ the terms of order 1 in A from those of order at least 
A; the former are very exphcit and can be put in a form similar to (even 
thought a bit more cumbersome than) the second line of Eg. (11. 151) : if we 
isolated them from R^"'\ then the remaining contributions to R^'^^ would be 
bounded by the r.h.s. of Eg. (II .16]) times |A|. Finally, in the massive case, 
m* 7^ 0, Eg. (ll.l6p can obviously be improved, by taking into account the 
exponential decay of the propagators on distances larger than (m*)^^. In 
this respect, the bound Eg. (ll.l6|) is stated with the massless case in mind, 
this being the most difficult to treat. 

1.2 Ising models, fermionic systems and a comparison 
with existing literature 

In order to put our results in the right context and clarify the relation be- 
tween spin and fermionic systems, we recall here some basic facts about the 
Ising model and its generalizations. Of course the literature on the Ising 
model is immense and we will give here only a partial view on the aspects 
more related to our work, pointing out the connections with the more recent 
literature and the main open problems. 

Let us start by discussing the relation between the integrable Ising model 
and fermionic systems. Schultz, Mattis and Lieb [52] derived the exact 
solution of the n.n. Ising model by a method alternative to the Onsager's 
one: they showed that the transfer matrix can be written as the exponential 
of a guantum Hamiltonian describing non interacting fermions on a one- 
dimensional lattice, and they computed its trace by using second guanti- 
zation methods. The Schultz-Mattis-Lieb solution points out a connection 
between two apparently unrelated systems, namely 2D classical spin systems 
and guantum many body ID non relativistic fermions. The relation between 
the critical Ising theory and a relativistic fermionic QFT can then be realized 
in several ways. One is to start from the Schultz-Mattis-Lieb representation, 
as in [32l |5T], and then use the fact, discovered by Tomonaga [55], that a 
many body system of non relativistic fermions in d = 1 can be effectively 
described in terms of relativistic fermions in 1 + 1 dimensions. Another pos- 
sibility is to start from the Grassmann integral representation of the Ising 
model partition function, as proposed by Hurst and Green [31j and Samuel 
[50] , see also Itzykson-Drouffe [32]. This approach starts from the obser- 
vation, see [13], that the generating function of the n.n. Ising model with 
periodic boundary conditions can be written in terms of four Pfaffians, each 
of which can be represented as a suitable Grassmann integral. Recalling the 
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definition of Z{A), Eq.f ll.7p . the outcome is (see e.g. [23150]) 

1 V 



Z°(A) := Z(A] 



A=0 



2^ TcZ'^iA) 

ae{±}2 



:i.l7) 



with Ti 



-T4 



1 and 

2 



5($,A) 



-2a 



Y\ ]^cosh(/3J + aAxj 
xeAj i=i 



a 1^ tanh(/3 J + aAx j)Ex^j + 



Here ifx, H^, Vy^, \4 are independent Grassmann variables (see e.g. for 
the basic properties of Grassmann variables and integration), four for each 



lattice site, and Ey^ i = Hy^H, 



x+aei ) 



while i?x2 



^x^^+ae2 ■ 



Moreover, 



$ = {-ffx, -f^x, V^x, Vy}y^(z\a^ denotes the collection of all of these Grass- 
mann symbols and P$ is a shorthand for Y\y:dH y^dH^dV y^dV^x.- The label 
ai,a2), with ai,a2 G {±}, refers to the boundary conditions, which 



are periodic or antiperiodic in the horizontal (resp. vertical) direction, de- 
pending on whether ai (resp. 0^2) is equal to + or — : 



UiHy 



^x+Le, '- 



= OiiHx 

aiVx . 



;i.i9) 



These boundary conditions may be important or not in the thermodynamic 
limit, depending on whether we are at or outside the critical point. If /9 7^ f3c, 
then |Z° (0)/Z° _(0)| converges exponentially to 1 as L — )■ 00, see [HI Ap- 
pendix G]; moreover, if /3 > f3c, then Z^{0) is positive as L — )■ 00 for all 
a, while if /3 < /3c, then Z^ _^{0) is definitely positive and the other parti- 
tion functions with different boundary conditions are negative. At /3 = (3^, 
Zl _^_{0) is exactly zero while the other partition functions can be expressed 
in terms of special functions [161 113] • 

The Grassmann representation is particularly convenient for computing 
energy-energy correlation functions: for instance 



'XJl 



T I 



,2\2 



2Z%0) 



{Ex J] Eyjl) 



ct,L ' 



:i.20) 



where {■)aL^^ average with respect to the Grassmann "measure" V^e^^"^'^^ 
with a boundary conditions; we shall also indicate by (■)^ the L ^ 00 limit 
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of mentioned above, if /3 7^ /3c then {E^^j; Ey')'^ is exponentially 

insensitive to boundary conditions as L — 00, so that, for all /3 7^ /3c, 

lim I = (1 - . (1.21) 

A similar formula is true for the m-point functions, m > 2: 

^im (e.,,,; ■ ■ ■ ; e^^,Jl^\^^^ = (1 - t')"^(i5;x„,,; ■ ■ ■ ; i?x„;,.)!,_ , (1.22) 

valid for all /3 7^ /3c. As mentioned above, if A = the Grassmann "mea- 
sure" used to compute the r.h.s. of this equation is gaussian (i.e. S'($, 0) 
is quadratic): therefore, the r.h.s. of Eq. fll.22p can be computed via the 
fermionic Wick rule, leading to a sum over all the Feynman graphs obtained 
by pairing ( "contracting" ) the Grassmann fields involved; truncation means 
that only connected Feynman diagrams should be considered. From a graph- 
ical point of view, Eq. fll.22p can be graphically interpreted as a sum over 
simple loop graphs, as in Fig{TJ The fermionic representation outlined above 
also allows one to compute the spin-spin correlations, even though these are 
represented by much more involved expressions. In particular, if x > 0, 

(a(.,o)ao,o>,,J,^o = t^((e»(*-^-*)&"^(^^o),.))^ , (1.23) 

where = J2a ^2z" (')q l- Eq.f ll.23p can be equivalently rewritten as 

(1.24) 

; -E'(z™,o),i))l • 

In other words, the spin-spin correlation can be expressed as a series of trun- 
cated energy correlations. Such representation is not very manageable and 
the asymptotic behavior of {o'(x,o)<^o,o) ij\x-o '^^ close to the critical point is 
usually derived by the (equivalent) representation in term of a determinant 
of an a; X a; matrix [13], even though several attempts have been pursued 
along the years to resum such series (see below). 

In any case, the above Grassmann representation looks quite different 
from those describing relativistic QFT fermionic models. However, we can 
perform a suitable change of variables [32] from {Hy^, i^x, Vy_, V^} to the crit- 
ical modes {V^x.it; Xx,±}, see next section for a precise definition; the prop- 
agator of the ip field is massless at the critical point, while the x fisld is 



+\m 2;-a 
m>l ^1=0 Zm=0 
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uniformly massive. After the integration of the massive modes, which can be 
performed exactly, we find that if A = 



(V'x,,.,---^x„,.J„,i = ^ / [ n #xJe^('^Vx,.,---^x„,.„ (1.25) 



xeAJ, 

U} = ± 



where Ma is a normalization and, if ipy, 



m --i I ^^ii {!:::% l^li)^. (1.28) 

where dj is a lattice approximation of the derivative in the j-th coordinate 
direction, A^, is a second order lattice differential operator, formally vanishing 
as a — 7- 0, and a = f *7*(o)°'' mass, where tc{0) = V2 — 1. Remarkably, 

the above functional integral coincides with the lattice regularization of a 
QFT describing Majorana fermions in ci = 1 + 1 dimensions, with a Wilson 
term (the second order operator A^) which allows one to avoid the fermion 
doubling problem [52] . 

Once again, the r.h.s. of Eq.f ll.25]) can be computed exactly in terms 
of the fermionic Wick rule, with propagator given by the covariance of the 
quadratic form in Eq. fll.26p . In the scaling limit, this propagator tends to 
0° of Eq. fll.9p : correspondingly, the r.h.s. of Eq. fll.22p . after the change of 
variables from {i^x, -f^x, ^x, K:} to {V^x,±,Xx,±} and the integration of the 
X fields, tends to the r.h.s. of Eq.f ll.Sp . If m*=0, the limiting theory is a 
scale-free relativistic Conformal Field Theory (CFT). Conformal invariance 
appears to be a robust property, stable under changes in the shape of the 
box and of the underlying lattice, as proved by Chelkak and Smirnov [15] by 
making use of suitable fermionic operators on the lattice [IB] and the ideas 
of Schramm-Lowner Evolution (SLE) [55] . 

While the discussion of the continuum limit of the energy correlations 
starting from the exact solution (that is, from the Grassmann integral rep- 
resentation) is straightforward and fully rigorous, the analogous discussion 
for the spin-spin correlation is much harder. One tempting route to its com- 
putation is to start from Eq. fll.24p and then replace the sums by integrals 
and the energy truncated expectations by their scaling limit Eq.f ll.Sp : this is 
the strategy followed by [21 |T71 158] . However, these replacements introduce 
spurious ultraviolet divergences that require a suitable interpretation; some 
justification of this procedure has been provided by Dotsenko and Dotsenko 
[T7] , but it is fair to say that so far there is no mathematically sound way to 
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resum the series in Eq. fll.24p to get the spin-spin correlation. A better alter- 
native approach is based on bosonization, which connects this (among many 
other quantites) to observables of a discrete height model [T9]. This height 
model has long been known to converge in a weak sense to the Gaussian free 
field [37j , and a recent preprint by Dubedat [18j extends this convergence to 
the relevant class of observables. 

Let us return now to the non integrable model Eq.f ll.ip with A 7^ 0. A 
basic principle in statistical physics, which has a wide experimental confirma- 
tion, is universality [51 [201 EH] , which tells us in particular that the physical 
properties of a system close to a second order phase transitions are largely 
independent of the microscopic details of the interaction among its elemen- 
tary components: only a few general properties related to dimensionality and 
the symmetries of the Hamiltonian matter as far as the computation of its 
critical exponents is concerned. It is widely believed that the Ising model 
with finite range interactions Eq. dl.ip belongs to the same class of universal- 
ity as the n.n. Ising model, that is, the critical exponents of the two models 
are the same. On the other hand, other physical quantities, like the critical 
temperature, are expected to be non-universal. However, a mathematical 
proof of this conjecture is very difficult: universality in two dimensions is not 
a trivial issue at all. A striking illustration of this fact is provided by the 
exact solution of the eight vertex model by Baxter [1], which came after the 
solution of the six vertex model by Lieb [lU], and showed that the critical 
exponents of these vertex models are continuous non trivial functions of the 
coupling. Eight vertex models can be represented equivalently as a pair of 
Ising models coupled by a four spin interaction [5]; therefore, the Ising uni- 
versality class should be stable under in-layer perturbations, but not against 
perturbations coupling two different layers. The classification of the possible 
critical theories close to the Ising model is by itself a very rich and inter- 
esting research field. It is based on the remark, due to Belavin, Polyakov 
and Zamolodchikov |6], that if a 2D critical theory admits a scaling limit, 
this should be invariant under the infinite dimensional group of conformal 
transformations of the complex plane. Such a large symmetry group imposes 
infinitely many constraints on the correlation functions and in some cases 
these are sufficient to compute them in closed form. Even more strikingly, 
in many cases the critical exponents are all functions of a single parameter, 
the central charge c, which in turn can be computed from the self-correlation 
of the energy-momentum tensor T{z), {T{z)T{0)) = ^ [HIES]- In our case, 
T{z) = —^ipdip and {T{z)T{0)) = from which c = | and the simplest 
conformal field theory with central charge 1/2 has two primary fields with 
critical exponents 2 and 1/4, respectively, which are naturally identified with 
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the energy and the spin operators of the continuum hmit of the critical Ising 
model [HI 123] ■ Universality and conformal invariance of the nearest neighbor 
Ising model on regular lattices and domains of different shapes has been re- 
cently proved by Smirnov [53], Chelkak and Smirnov ^] and by Chelkak, 
Hongler and Izyurov [H], by SLE methods. However, it should be stressed 
that the results of Smirnov, Chelkak and collaborators heavily rely on the 
underlying integrability properties of the nearest neighbor Ising model and 
are very fragile under perturbations of the form Eg. (11. II) . The problem of 
proving the existence of the scaling limit of the critical theory Eg. (11. II) at 
A 7^ and of its conformal invariance still remains a big challenge. 

A different approach to universality which is suitable to include pertur- 
bations of the form Eg. Ol.ip is provided by the Renormalization Group [56] . 
A rigorous application of this idea to the present context has been proposed 
by Pinson and Spencer [171[5lj. As mentioned above, their starting point is 
an exact expression of the partition function of certain non-integrable Ising 
models in terms of a non-gaussian Grassmann integral, which can be stud- 
ied by constructive QFT methods, similar to those used by Lesniewski jSH] 
to analyze the Yukawa2 QFT: in both cases the interaction is irrelevant in 
the Renormalization Group sense. Using such methods, Pinson and Spencer 
proved that the critical exponent of the energy correlation is independent of 
the next-to-nearest-neighbor interaction, provided this is chosen of a suitable 
form. Their proof provides the first example of universality of a critical ex- 
ponent in a perturbed 2D Ising model. This approach based on fermionic 
mapping and constructive Renormalization Group analysis can be applied to 
coupled Ising models, like the Eight Vertex or the Ashkin- Teller model; in 
such cases, as shown by Mastropietro [H], the interaction becomes marginally 
relevant in the Renormalization Group sense and the energy exponents are 
non universal continuous function of the coupling. Extensions of such meth- 
ods allowed Benfatto, Falco and Mastropietro to prove the Kadanoff relations 
[3^ between the specific heat, the energy and the crossover critical exponents 
in eight vertex and Ashkin- Teller models |8]; they also allowed Giuliani and 
Mastropietro to compute a new critical exponent controlling the crossover 
from universal to non-universal behavior in the asymmetric Askhin- Teller 
model |27j . It should be remarked that these exponents cannot be computed 
by other means, since the Ashkin- Teller model is not solvable (and we also 
recall that the eight vertex model is solvable but only certain exponents can 
be deduced from the solution). 

In this paper we extend the analysis of Pinson and Spencer [171 in 
several directions. First, we generalize the class of spin perturbations that 
can be considered: while they reguired special next-to-nearest neighbor in- 
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teractions, we just need the function v in Eq. fll.ll) to be of finite range and 
symmetric under the natural lattice symmetries; the resulting exact fermionic 
action is defined in terms of an exponentially decaying self-interaction of the 
involved Grassmann field, see Proposition [1] in Section 12.11 below. It should 
be noted that even though the interaction in terms of spins is finite range, 
the corresponding fermionic interaction is of infinite range and of arbitrary 
high degree in the fermionic field, but exponentially decaying on the scale 
of the lattice. The result is based on a first cluster expansion, which is ex- 
ponentially convergent provided that the strength A of the perturbation is 
small enough. Second, we combine the analysis of Pinson and Spencer with 
one used by Benfatto, Falco and Mastropietro to prove and control the 
convergence of its energy correlations to the continuum limit. In this way, 
besides the critical exponent of the 2-point energy correlation, we can ex- 
plicitly compute the scaling limit of all the multipoint correlation functions: 
in particular we have a constructive procedure to compute the amplitude of 
the energy correlations and the subdominant corrections that vanish in the 
scaling limit. As a technical point, let us also mention that, as compared to 
[7], our bound on the subdominant corrections is optimal from a combinato- 
rial point of view: it grows as m! as m — )■ oo, exactly as the dominant term, 
contrary to the bound in [7], which grows as (m!)", a > 1. Our results can 
be seen as a strong statement of universality with respect to perturbations of 
the form Eq. fll.ip : the critical correlations in the scaling limit have the same 
analytical expression as those of the n.n. Ising model, up to a renormaliza- 
tion of the wave function and of the mass. 

Of course, many important problems remain to be faced. One is to repeat 
the analysis for the energy correlations on the torus, taking the scaling limit 
directly at the critical point. However, the most urgent problem is to prove 
universality of the spin critical exponent. As we have seen before, the explicit 
expressions for the energy correlations of the Ising model have been used by 
Dotsenko and Dotsenko [TTj to compute the spin-spin critical exponent (mod- 
ulo a number of audacious exchange of limits). The spin-spin correlation for 
the model Eq. (11. ip with A 7^ can be still expressed by a formula similar to 
Eq. (ll.24p . with {t~^ —t) replaced by a renormalized coefficient tx and the free 
truncated energy correlations replaced by the interacting ones at A 7^ 0. If we 
proceed as in [T7] and, in the scaling hmit, we replace the sums in Eq. (ll.24p 
by integrals and the truncated energy correlations by their limiting value, we 
get an expression essentially identical to the non interacting one, with the 
important difference that the m-th order term in the series has a prefactor 
proportional to {txZ{X))'^ . After a resummation (if formal) of the series, the 
factor t\Z{\) appears to be related to the critical exponent of the spin-spin 



16 



correlation. In this perspective, the proof of universahty for this exponent 
could be reduced to the proof that the combination txZ{X) is identically 

1 as a function of A. However, one might worry that the exchange of the 
scaling limit with the integrals could produce an extra finite renormalization 
that could be difficult to control; this question requires further investigation. 
Another possible strategy for determining the exponent of the spin-spin cor- 
relation is by computing the self-correlation of the energy-momentum tensor: 
if we could identify the right lattice counterpart of this operator (which is 
usually defined directly in the scaling limit) then we could have direct access 
to the central charge of the interacting theory (which should be 1/2 if the 
universality principle is correct). 

The rest of the paper is devoted to the proofs of Theorems 11.11 and II. 2[ 
In Section [2] we prove the Grassmann representation of the generating func- 
tion for the multi-point energy correlations. In Section [3] we describe the 
Renormalization Group procedure used to control the Grassmann generat- 
ing function uniformly in the mass m*. Finally, in Section H] we explain how 
to adapt the general expansion and the bounds discussed in Section [3] to the 
multi-point energy correlation functions and conclude the proof of the two 
theorems stated above. 

2 The representation of the interacting Ising 
model in Grassmann variables 

2.1 The generating function for the energy correla- 
tions 

Consider the model Eq. fll.ll) where f(x) has the properties spelled after 
Eq. fll.ip . Without loss of generality we can assume that the additional near- 
est neighbor interaction is zero: v{aej) = 0, j = 1,2. The interaction is 
written in terms of the macroscopic coordinates for a unified notation, but 
in fact depends on the lattice distance; that is f (x — y) = Vo{^^-^) for some 
vq independent of a. We shall assume that v is normalized in such a way 
that 1^(^)1 ~ -^^^ notational simplicity, we shall also assume that 

a = io2~^ , where io is the macroscopic unit length and N & N. Moreover, 
we define B^j to be the set of nearest neighbor bonds in = ci'^'m- 

As reviewed in the previous section, in the simple Ising model case (A = 
0), the generating function for the energy correlations can be represented 
in terms of a gaussian Grassmann integral, see Eqs. fll.17p - fll.18p . A similar 
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(2.2) 



representation in terms of a non-gaussian Grassmann integral is valid also in 
the A 7^ stated in the following proposition and proved below. 

Proposition 1. There exists Ao > such that, if \X\ < Xq, then for any 
m-tuple of distinct pairs (xi,_7i), . . . , (x^, j^), with m > 2, Xj e and 

e{l,2}, 

am 

where S( A) is the Grassmann generating functional 

where t :— ta,nh{/3J), {E, A) :— Yl^=i J c^^-^xj^x,j and: 

• Cm is a normalization constant, defined as 

Cm = {-2a-' cosh\pj))^\''^^^^ J] cosh2(^^;(x - y)) (2.3) 

{x,y} 

with Vm(A) an analytic function of \, independent of a and satisfying 
the hound \ Vm{X)\ < C|A|M^, for a suitable C > 0; 

• St{^) is the unperturbed quadratic part of the action, defined as 

St{^) = aJ2{tE^,i+tE^,2+ (2-4) 

• V($, A) is a polynomial in {-ExjjxeA" '^^^ {^xjIxga^ ? which can be 
expressed as 

V($, A) = J] 5Z / • • • • • • ^^"^ • (2.5) 

n,m>0 ji,---,jn 
n+m>l 



•M^jj'(xi, ■ ■ ■ , x„; yi, . . . , y^) JJ Ex.j, H ^yi"^-' 



i=l i'=l 



w^ere j = (ji, . . . , j„), / = . . . J'^) and, z/ x = (xi, . . . , x„) and 

y = (yi,---,ym); 



g K5(x,y)/a 



\Wjji^;y)\ < c-+n^(^|A|)-^^{^^^(^+H}a-(2-n-m) (2.6) 



Q2(n+m— 1) 
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for suitable constants C,c, n > depending only on Mq ( the range 
of the interaction); here 5(xi,...,ym) is the tree distance of the set 
X = {xi, . . . , Ym}, that is the length of the shortest tree graph composed 
of bonds in B%j which connects all the elements of X. 

Remark. The factor ^2(n+"l^i) in Eq. fl2.6p is normalized in such a way 
that its Li norm is essentially independent of a, that is 

I r g-K<5(x,y)/a 

where is a constant independent of a. On the other hand, the factor 
(/3|A|)™^^''^^''^^"'"'"'"^^ measures the smallness in A of the kernel W ^ while the 
factor a"*^^"""™) (or, better, its a-dimensional version, (a/£o)~^^~"~™^) plays 
the role of its scaling dimension. Recall that a/^o = 2~^, see the beginning 
of Section [2m therefore, the scaling dimension of the kernel can be rewritten 
as 2^(^~"~™') . We shall see in the following that the generating function can 
be integrated by an iterative multiscale procedure; at each step, it will be 
expressed in terms of an effective potential on scale /i, with h < N, analogous 
to V($, A), whose kernels W^'^^ have decay properties analogous to W, with 
the important difference that the scale 2^ is replaced by 2^ and, therefore, 
the scaling dimension 2^(2""""^) is replaced by 2'^(^~"'~™'). The result of the 
iteration, as /i — )■ — oo, gives the generating function of interest. The relevant 
scaling properties of the multi-point energy correlation function as a — )■ and 
/3 — > /3c will be controlled in terms of the kernels W^'^'' and of their limits as 
— 7- oo and h — )■ — oo. 

Proof of PropositionUi For notational convenience, given a bond b G Bij, 
we denote by e^, Ei, and Aj, the corresponding bond operators: that is, if 
b = (x, X + aei) (resp. b = (x, x + 062)), then Eb = £x,i, Eb = -Ex,i and 
Ab = Ax,i (resp. Eb = £x,2, Eb = -E'x,2 and Ab = v4x,2)- The key to the first 
step in the proof is the remark is that if 61, . . . , 6„ are n distinct bonds, then 
the Grassmann representation for the nearest neighbor Ising model induces 
the following: 

{-2a-')''" [ n cosh(/3 J + aAb)] ■ (2.7) 
J D$ {a-Hb, + (1 - tl)Eb,) ■ ■ ■ {a-Hb„ + (1 - t^jE,Je^(*'^) , 
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where tf, = tanh(/3J + aAh). This correspondence is invahd for repeated 
bond variables: note that [a~Hh + (1 - tl)Eb]'^ = a~Hl + 2a-^{l - tD^Eh, 
while el = a~^. This last observation can be used to remove repeated bond 
operators from any expression; therefore, in order to derive a Grassmann rep- 
resentation for Z{A), it is enough to express the interaction term (i.e. the 
A-dependent term) in Eq. fll.7p as sum of products of distinct bond operators; 
then we can replace every bond operator Eb by a~^th + {l — tl)Ef,, in the sense 
explained above, and finally we can re-exponentiate the big sum of products 
of Grassmann variables, so obtaining the desired Grassmann functional inte- 
gral representation of the Ising model at hand. This can be implemented as 
follows. By definition, 

ctGHm {x,y} 

Consider a pair of sites {x, y} contributing to the product in the r.h.s. of 
this equation, i.e., a pair of sites such that |x — y| < Rq. Note that (TxCTy 
can be rewritten in terms of a product of energy density operators localized 
along a path connecting x and y on the lattice: o"xO"y = i[/x,y + |-Dx,y, where 
U^,y = tlb€Cu{^,y) and Dx.y = UbeCn(^,y) ^^b- Here Ct/(x, y) and Cd{^, y) 
(where U and D stand for "up" and "down") are the two paths connecting 
x and y on the lattice described in Fig. [31 Note that we choose paths in this 
way in order to be sure that we have an expression which manifestly retains 
the rotation and reflection symmetries of the original interaction. In terms 
of these "string operators" , we can use the identity 

e^^ = coshx(l ± tanhx) 

and the fact that U,D = ±1 to rewrite the product in the r.h.s. of Eq. fl2.8p 
as: 

JJ g/3A...(x-y)., ^ -Q ^M.(x-y)(f/.„+D.„) ^ J -Q cosh2(f y(x-y))' 
{x,y} {x,y} {x,y} 

[II (l + tanh(i/3A^;(x - y))U^,y) (l + tanh(i/3A^;(x - y))/^x,y) ■ 
{x,y} 

The second line is a product of binomials, each consisting of 1 plus a non 
trivial term; these non trivial terms can each be graphically associated with 
a "string" S (i.e., the union of the bonds in C[/(x, y) or Cj:)(x, y)), either 
of type U or D, depending on whether it is associated to the path Cf/(x, y) 
or Cd(x, y); in both cases we shall write vs '■= v(x — y). Note that the 
assumption that the interaction has range oMq means that all of the strings 
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(b) 

Figure 3: The two paths C[/(x, y) (on the left) and C£)(x, y) (on the right) 
in the two cases where: (a) xi < yi and X2 < y2', (b) xi < yi and X2 > y2- 
If Xi — yi or X2 — y2: we choose C[7(x, y) = Cd(x, y) to be the horizontal or 
vertical segment connecting x with y; this means that horizontal or vertical 
strings are counted twice. 



we are considering consist of no more than Mq bonds. If we now expand the 
product 

[ II (l+tanh(i/3A^;(x-y))C/x,y)(l+tanh(i^AT;(x-y))D,,y) 

we get a sum over all subsets of the collection of strings which appear in the 
interaction under consideration, in which each term is a product of string 
operators. Moreover, every such subset can be thought of as a collection of 
its maximal connected components (here we say that a set of strings S = 
{Si, ... , Sm} is connected if, given 1 < io < Jo ^ we can find a sequence 
(S'io, Sij^,. . . ,Sip = Sj^) such that n Si^_^^ ^ 0). Prom a graphical point of 
view, every connected component S corresponds in a non-unique way to a 
polymer 7(5), i.e., a connected set of bonds. It is helpful to color the bonds 
in 7(5) black or gray, depending on whether the given bond belongs to an 
odd or even number of strings in <S, and denote the set of bonds thus colored 
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black by bl(iS). We call the collection of terms associated with a polymer 7 
the activity of 7 and denote it by 2; (7) and noting that a^ef = 1 

^(7)= Yl [lltanhC^pXvs)][ H aet] , (2.9) 

5 connected: 5g5 feebl(<S) 
7(5)=7 

in terms of which we can finally rewrite 

-Q g/3Aa..(x-yK = J Y[ cosh^ (i/3At;(x - y))] 5^ V'(r) Y[ 2(7) (2.10) 

{x,y} {x,y} rcA 7Gr 

where the sum ^pcA r.h.s. runs over sets of polymers, T = {71, . . . , 7„}, 

to be called polymer collections, such that each polymer is contained in A = 
A^, i.e. it is formed by bonds in B^j. Moreover, the function ipdji, . . . , 7n}) 
implements the hard core condition, that is (p is equal to 1 if none of the poly- 
mers overlap, and otherwise; the term with F = should be interpreted as 
1. Note that the r.h.s. of Eq. f l2J0|) . like that of (E9]), is multilinear in the 
bond variables (i.e., each bond variable Sb appears at most once in every term 
in the sum). Hence, plugging Eq. fl2.10p into Eq. fl2.8p and using Eq. fl2.7p . we 
find: 



ZiA) = J]^(-2a-2)^"'[ n coshiP J+aAb)] [ J] cosh' {1/3 Xv{^-y)j 

hdBli {x,y} 

■ f V^Y.^{V)[\[Ul)Y^'''''\ (2.11) 



rcA 7gr 



where 



Cg(7)= W^^^KlP^vs) n {h + a{l-tl)Eb) , (2.12) 

S connected: S^S feebl(5) 
7(5)=7 

This is exactly the sort of expression which is the subject of the standard 
cluster expansion (for a presentation in a convenient form, see [221 Chap. 7]), 
which gives the identity 

E ^(r) n Cg(7) = exp { 5^ v^^(r) n Cg(7)} (2.13) 
rcA 7er rcA 7er 

which is valid provided that the sum in the r.h.s. is absolutely convergent. 
In the r.h.s. the sum over F involves polymer collections that include over- 
lapping and even repeated polymers (we let F(7) be the multiplicity of 7 in 
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— )• tanh(|/3At'5j tanh(|/3At'52) tanh(|/3Af53, 
{K + (1 - t K + (1 - tl)aE,,) (t,3 + (1 - t^3)aE,3) (t,, + (1 - tlJaE,,) 

Figure 4: Example of a set of strings S = {5*1, 5*2, 5*3} and the corresponding 
coloring of the associated polymer 7. The contribution to Cg associated with 
S is given below. 



r) and y9-^({7i, . . . ,7„}) are Mayer's coefficients, which admit the following 
explicit representation. Given 71, . . . , 7^, consider the graph Q with n nodes, 
labelled by 1, . . . , n, with edges connecting all pairs i, j such that 7^ fl 7^ 7^ 
{Q is sometimes called the connectivity graph of the collection of polymers 
7i, . . . , 7„). Then one has = 0, V^"^({7}) = 1 and, for n > 1: 

^^({71, . . . , 7n}) = ^ ^(-l)number of edges in (2.14) 

■ ccg 

where F! = 117^(7)! and the sum runs over all the connected subgraphs 
C of G that visit all the n points l,...,n. In particular, if n > 1, then 
V9^({7i, . . . , 7„}) = unless {71, . . . , 7„} is connected. 
Plugging Eq. fl233|) into Eq.f l2TT|l gives 

Z{A) = ^^(-2a-2)A/^J -Q cogh(/3J + aA)] • (2.15) 
■[ W cosh2(i/3At;(x-y))] j p$e^(W^(*-^) , 



{x,y} 



with 



V($,A):=^y,^(F)nCG(7), (2.16) 

rcA 7er 

provided that the r.h.s. of this equation is absolutely convergent, in the fol- 
lowing sense: V($, A) is a polynomial in the Grassmann variables {-ffx, H, 
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Figure 5: A fully colored contour as used in evaluating W: bi is orange, 62 
is red, 63 is yellow and 64 is black; the other unlabeled bonds are grey (color 
online). By definition, the activity of this colored contour is proportional to 
(-2t(l - t^)) aEb,aAb^ x (1 - t'^)aEb, x (1 - t^)aAb^ x t. 



Vx, Vx}, whose coefficients can be written as infinite sums over polymer con- 
figurations {71, . . . ,7„}; by absolute convergence of V, we mean that each 
of these infinite sums are absolutely convergent. As far as the dependence 
of V on A is concerned, we note that since Ceil) is defined through finite 
sums and products of hyperbolic tangents U = tanh(/3J-|-ay4;,), the activities 
Ceil) are analytic as functions of each Ab. However, only the terms at most 
linear in each Ab can contribute to the computation of (eb^; ■ ■ ■ ; £bm)'^ with 
61, ... ,6m all distinct. Therefore, to the purpose of computing these trun- 
cated expectations, we can safely replace Cg(7) by its multilinear part in 
{Ab}bGBij, namely: 

Cg(7) ^ Coil) := E E ^' ^) n n "^'^ ' (2.17) 

RCj YC-y beR beY 

with 

ft, b^RUY 
C(i?,F;7):= E n*^^^(i'^^^^) n { 1-^'' be RAY 

S connected: S£S feGbl(5) I —2t(l — t^), b E RHY 

7(5)=7 
h\{S)2RUY 

(2.18) 

where {R, Y; 7) can be associated with a decorated contour where the bonds 
in R are drawn red, all bonds in Y are drawn yellow, and all bonds in 
both are drawn orange (see Figure Similarly, we can replace the fac- 
tor ribeH" cosh(/3J + aAb) e*^*^*'"*^ in Eq. (l2.15p by its multilinear part in 
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{Aiy}h<zi3ij. Performing these replacements corresponds to replacing the gen- 
erating function in Eg. (12 .15]) by: 

Z{A) ^ e'^'^^^'EiA) , (2.19) 

where S(A) is the same as in Eq. (l2.2p . with 

Vm{X) + v{^,a) = ^^^(r)nCG(7) = 

rcA 7Gr 

= E E ^) [ n ^^'-l [ n ' (2-20) 

and Vm(A) = 1^(0,0). Here the coefficients W{R,Y) are defined as 
W{R, y)= J2 ^ ^^(r(7)) Yl n 7.) , (2.21) 

">0 ■ (iJi,...,_R„)GPfl,^ j=l 

7=(7i,...,7„) (Yi,...,y„)G-Py,^ 

where 7 = (71, . . . ,7„) is an ordered n-tuple of polymers, r(7) is the cor- 
responding unordered n-tuple and the combinatorial factor is used to 
pass from the summation over F to the one over 7. Moreover, Va^-f is the 
set of ordered n-ples of disjoint sets (-Ri, . . . , i?„), such that Ri C 7^ and 
U^^i-Rj = R- Once again, these rewritings are valid provided that the sums 
entering the definition of W{R, Y) are absolutely convergent. Note also that 
the prefactor e"'*^>>^b i^j^g r.h.s. of Eq. (l2.19p is irrelevant to the purpose of 
computing the truncated expectations (efe^; ■ ■ ■ ]£b„i)'pL with m > 2, which 
explains why we did not insert this prefactor in the statement of Proposition 

m 

The absolute convergence of Eg. (12.2 ip follows from an important result 
in the cluster expansion [221 Proposition 7.1.1], which states (in part) that 
for any positive function on polymers which decays as 

fil) < z/o'l^le-2''°^(^)/" (2.22) 

for some kq > and uq sufficiently small (here I7I is the number of bonds in 
7), then 

sup Yl ^''(r)! n/W^2i.oe-^^/". (2.23) 
r={7i,...,7„} 76r 

diam(Uj7j)>-R 

Actually, the proof of [221 Proposition 7.1.1] implies a slightly more general 
bound; namely, given B C B%j, 

E n/(7) ^ C'i^o"""e-'=^''°'^^)/'^ , (2.24) 

r={7i,...,7„} 7Gr 
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for two suitable constants Ci,Ci > 0. 

To apply this, we first note that the generic contribution to the r.h.s. of 
Eg. (12.211) vanishes unless the collection 7 = (71,..., 7„) is connected and 
touches all the bonds in B := RUY. Therefore, 



\W{R,Y)\< Yl ^lv'^({7i,---,7n})| Yl niC(^.'^.S7.)l 

(2.25) 



n! 

n>0 (i?i,...,i?,„)ePji,^ j=i 

■7=(7i,...,7„) (Yi,...,Yn)eVY,-y 



where 

n n n 

(Rl,...,i?„)GPfl,^ i = l i = l RjC-yj j = l 

(Yi,...,y„)GPy,-, Yj^-rj 
and 

/(7):=4l^lmax|C(/?,F;7)l (2-26) 

RC-y 
YC-y 

plays the role of the function /(7) in the l.h.s. of Eq. (l2.24p . In order to 
use the bound Eg. (12.241) . we need to prove Eg. (12. 22 1) . Using the definition 
Eg. (12181) of C(^, 7), it is clear that 



/(7) < Yl Hi/^i^^^i 

S connected: 5G<S 
7(5)=7 

m>iM be7 S9b 

— aMQ 

where the * on the sum Yl*SBb indicates the constraint that b is either 
the first or the last bond in S. By using the normalization | ~ -'- 

(which means that | Y^s^b l^sl = 1) and defining mo = mo(7) = ^ = JJ, 
we get: 

m>mQ 

<- ^"S(^)'"T^(^)" (-) 



mn! \ I J — ' m! 

m>0 

< (4e)l7l(^M)"^°e|/3|A||7l ^ ^2|7|^-2.o^(7)/a 
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which proves Eq.([223) with vl = 46^+5/31^1 (M)2iOo and e'^'^o = (^) ^'^^o . 
This concludes the proof of the representation Eq. fl2.2p . The bound Eq. fl2.6p 
on the decay of the kernels of V($, A) is just a restatement of the decay bound 
on W{R, Y) that we just derived (note that the factor a"*^^"""™-* comes from 
the factors a in front of the fields Eh and Aj, in Eq.f l2.20|) and from the 
definition of / dx, which brings along a factor a^; in fact, J rfx = XlxgAj^)- 
All that remains is to prove the stated properties of VmW = H^(0,0), 
where 



w^(0,0) = 5^¥.^(r)nC(0,0;7) 



(2.29) 



rcA 



7Gr 



supp r I 



nc(0,0;7) 



supp TBb 



7er 



where suppF = U^gr7 and | suppFj is the number of bonds in suppF. The 
sum in square brackets is independent of b, by translation invariance. More- 
over, it is independent of a, as it follows by the definition of C(0,0;7) and 
by a relabehng of the lattice spacing a. Therefore, W{^, 0) can be bounded 
by 2M^s(A), where s(A) is a bound on the expression in square brackets. A 
repetition of the argument used for W{R, Y) above implies that s(A) is of 
order A, which concludes the proof of Proposition [1] □ 



2.2 Majorana form of the action 

As we have shown, the non-integrable Ising model under consideration can 
be expressed in the form of an interacting fermionic system, described by 
the action Eq.f l2.2p . which consists of a leading term, S't($) + [t + (1 — 
t'^)Eb\Ab, plus an interaction, which vanishes as A — )• and, in this respect, 
is "subdominant" . Since in the following we want to treat this subdominant 
term as a perturbation in the vicinity of the critical point, it is convenient 
to use coordinates adapted to the critical modes of the leading term, as 
described in the following. 

The Fourier transforms of the Grassmann variables are defined as: 

:= Yl ' := Yl ' (2.30) 

V^:=a^ Y e^''"^, K := Yl ^^'^'^x , (2.31) 
where, if a = (ai,a2), then k e V^j := {^(n + |a) : n G ZVMZ^}. The 
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inverse transformation reads 

-y 



— ik x TT 



7^ E e-^'^'^^k , (2.32) 



ke©? 



jkx 



(2.33) 



If •I'k is the column vector with components iJk? -f^k, ^k? ^k? respectively, the 
leading quadratic part of the action, St{^), can be rewritten in Fourier space 
as: 

S^i^) = 4 E "^"^ka^k , (2.34) 



ke©? 



where 



k •- 



v 





-1 - te*'"''! 
1 
1 



1 + te 


-1 
1 



-iaki 



-1 
1 



-1 -te'"''^ 



-1 \ 

-1 

1 + te^'''''^ 

/ 



(2.35) 



the inverse of which has the meaning of free propagator of the Grassmann 
field $. Note that C^^ is singular only at k = and t = \/2 — 1 (criticality 
condition). In this case, Cq has two vanishing eigenvalues and two purely 
imaginary eigenvalues ±ia^^\/2, and the corresponding eigenmodes read: 



u 





/ e^f 


e-*f 1 


-A 


1 


-j- 
e 4 


1 


z 


2 


j- 

—6 4 


-e-f 1 


—i 






-e^f 1 


^1 



(2.36) 



The natural variables at the critical point are the "critical eigenmodes" de- 
fined by the unitary transformation U in Eg. (12. 36 p . namely 



Xk,+ 



(2.37) 



For later convenience, we rescale these variables as (where w = ±) 

V'L = ^^k,^ , xlo; = -^Xk,<^ , (2.38) 
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so that, defining ipi^ (resp. Xk) as tlie column vector witli components 
tpk,+ ,^pk- (resp. Xk,+5Xk,-), we can rewrite: 

Sm = Yl (^-kC^^(k)4 + X-kC^(k)xk) + X) , (2.39) 

where, if ^ = ip, x- 



(a-^Hsmah + sinah) .a#(k) ^ \ 

^ Y —ta#[k) a [—tsmaki — smak2)J 

and 

(T^(k) = - ^cosa/ci + cosaA;2 — 2-^^^— ^ j , (2-41) 

a^(k) = - ^' cos afci + cos afca + 2 ^"*" . (2.42) 

Moreover, 

Q(^, X) = ^ E (^-kQ(k)Xk + X-kQ(k)4) , (2.43) 

n(h-\ — ( ^^^(~^sina/ci — sinaA;2) a~^(2cosa/ci — 'icosaA;2)\ 
ya^-'^(— icosafci + zcosaA;2) a^^(— zsina/ci + sinafc2)/ 

In terms of this notation, we can rewrite Eg. (12. 2 p as 

H(A) = ^A^-" / i^c.(#)Pc.(c^x)e^(^'^)+^('^'^'^)+^('^'^'^) , (2.44) 

where: 

• The normalization constant CM,a is defined as 

Cm,c. = ^(a^M)^*^V,Ar^,«Ar^,« , 

where: Cm was defined in Eq.([23D; the factor {a^Mf^^ takes into 
account the change of variable from the set of Grassmann variables 
{ifx, -f^x5 ^x} to {'?/'k,Lj? Xk,w}: in fact a computation shows that 

/P<f [ n ^ki^kFkFk] = a«*^^[detk,x(e^'^-)]' = (a^M)^^^^ ; 
kec" 

Tq, was defined after Eq. fll.l7p : M^^a and A/'^,a are the normalization 
constants of the two Grassmann gaussian integrations Paidtp) and 
Pa{dx), see next item. 
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The Grassmann gaussian integrations Pa{dtp), Pa{dx) are defined as 
PUd^lj) := ^ [ n n ^^k,.] exp { - ^ 5^ ^^i,C^(k)4} 



and J^-,p,ct,J^x,ct two normalization constants, fixed in such a way 
that / Po(#j = / Pc{dx) = 1. 

The source term B{^p, Xi ^) is the rewriting of (1 — t'^){E, A) in terms 
of the new variables, namely 

B(,^,X,A) = -, — J-. (2.45) 

^x,i(^x,+ - #x,- - Xx,+ + «Xx -)ri(-z^/'x,+ + ^x - + «Xx,+ - Xx -) 

+^x,2(^x,+ - V^x - + Xx,+ - Xx,-)r2(V'x,+ + V'x - + Xx,+ + Xx - 

where tj is the translation operator that shifts by one lattice step the 
argument of the field which it acts on: Tjtpy^^^ = V'x+aej.w and similarly 
for X- 

• V{ip, Xi ^) is the rewriting of V($, A) in terms of the new variables. 
It is easy to check that its kernels satisfy the same decay estimates as 
those of V($, A), see Eq. fl2.6p in Proposition [TJ 

As mentioned in the introduction, we are concerned with a scaling limit such 
that we take the thermodynamic limit M — )■ oo first, keeping /3 = /3(a) 7^ 
/3c(A), and then (simultaneously) a — )• and {3 — )■ /3c(A). In doing so, the 
resulting multi-point energy correlations are insensitive to the Grassmann 
boundary conditions, labeled by the four possible values of o:; this is true 
both in the A = and in the A 7^ case, see [41i Appendix G]. Therefore, 
the multi-point energy correlation functions in the specific scaling limit that 
we consider are the same as those computed from the following generating 
function: 

P(rf^/')P(dx)e'^^^'^)+^^'^'^'^^+^^^'^'^^ , (2.46) 

where P{dip) is a shorthand for P^^^{dip) and similarly for P{dx)- For future 
reference, let us note that the propagator of the ip and x fields associated 
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with the gaussian integrations P{dip) and P{dx) are given by, if -D^ (k) = 
a~^(zsina/ci ± sinafc2) and |£)a(k)p = —D^{k)D~{k), 

atA^-y) = j P{d^)i^.,.i^y,.' (2.47) 

27r 

and 



L^,. l^a(k)P + K(k)]2 ^-.a^(k) /^-(k) j , 



<.'(x-y) = j Pidx)x.,.Xy,.' (2.48) 

^ 27r >^ e-"-(-y) / D+(k) ta^{k)\ 

J^^^ + K(k)]2 {-^a^{k) D-{k)l^^, 

where Vm is a shorthand for T>^j~ . 

Remark. We define the unperturbed scaling limit as follows. Let /3(a) be 
fixed in such a way that, if t = t{a) = tanh (/3(a) J) and tc(0) = ^2 - 1, 

t{a) ~ 2 ' ^ 

where (j°(a) 7^ for all a 7^ and Ymva^o (t^ {a) = m*. The limit as a — ?• 
with t{a) fixed in this way will be referred to as the unperturbed scaling 
limit. The explicit form of the propagator of the ip field shows that in the 
unperturbed scaling limit it behaves as: 

\ 0/ N f dk e~'^^ fiki + k2 im* \ 
/(X) ^ B»(x) ^^^^ _ ^ j . (2.50) 

Note that the limiting propagator 0°(x) is normalized in such a way that in 
the massless case, m* = 0, it reduces to Eg. (11. lip . In the same limit, the 
propagator of the x field 5'^(x) tends to zero for every fixed x G M^; moreover, 
the combination f d^g'^{'x) tends to —c^a2, where = Imia^o o,a^{Q) = 

p|jy(tc(0) + V2 + 1) and (T2 = is the second Pauli matrix. In this 

sense, as a tends to zero: 

g^{x) ~ -aCy^5{x.)(T2 , (2.51) 
where 6{x.) is the Dirac delta function. 
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The quadratic coupling Q{iIj,x) iii Eg. (12. 46 p can be eliminated by the 
linear change of variables 

Xk^Xk + C-i(k)g(k)4, (2.52) 

which leaves S_ _(A) invariant. Note that for small k, the kernel C~^(k)(5(k) 
associated to this transformation is small, namely C~^(k)Q(k) = 0(a|k|); 
that is, from a dimensional point of view, the action of C^^(k)Q(k) on -i/^k 
is the same as the action of the differential operator ac^x- After the transfor- 
mation Eq.f l2.52|) we can rewrite: 

S_,_(A) = ^ / P(#)P(c/x)e^(^'>^'"')+^('^'^'^) , (2.53) 

where Bi^ip, x, A) and V{ip, x, A) are the rewritings of B{ip, x, A) and V{ip, x, A), 
respectively, in terms of the new variables; moreover, the gaussian integration 
P{dip) is defined as 

p(#) = ^[ n n ^^k,.] exp { - ^ v^-kC^^.(k)4} , (2.54) 



M 



where the normalization constant AT^ is chosen in such a way that / P{di)) = 
1 and 

C^k) = C4k) - Qik)C-\k)Q{k). (2.55) 

Remarks. 

1. It is easy to check that the change of variables Eq.f l2.52|) does not affect 
the bounds on the kernels of the effective potential; i.e., the kernels of 
B and V satisfy the same bounds as the kernels in Proposition [H see 
Eq.dlSD. 

2. The correction C^{k) — C^{k) is small at small k, i.e., 

Q{k)C-\k)Q{k) = 0{a\k\^) . 

In particular, this means that the unperturbed scaling limit of the 
propagator ^^^^/(x) = J P{dilj)ip:^^^ipo^^/, in the sense of the remark 
after Eq. 02.48 p . is the same as the one of ^''''(x), that is, it is equal to 
0°(x). 

3. The representation Eq. fl2.53p is valid also when A = 0, in which case 
V = and the scaling limit of the energy-energy correlations can be 
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computed explicitly. In fact, using the explicit forms of Bltp, x, A) in 
Eq. fl2.45p . of the linear transformation Eq.f l2.52]) and of the gaussian 
integrations P{dip)P{dx), we find that for all m-tuples of distinct points 
xi, . . . ,Xm, m > 2 the unperturbed scaling limit (in the sense of the 
remark after Eq. fl2.48p ) of the energy correlations is: 

lim lim (^xi i, ; ■ ■ ■ ; i )'^t \ t = 

= lim lim a-^'"— S_ _(A)L = (2.56) 



—I 



)m 
^g^(^xi, + ^xi,-; ■ ■ ■ ; V'x^, + V'x^,-) 



where in the last line £jo indicates the truncated expectation with re- 
spect to the gaussian fermionic integration Pgo{dip) with propagator g°; 
i.e., given n functions Xi, . . . , X„ of the Grassmann variables tp: 



El{X,- ■ ■ ■ = —log / Pgo((i^)e^^^^+-+^"^"l 



Qn 

(2.57) 

and a similar definition is valid for a more general gaussian fermionic in- 
tegration. In the last line of Eq.f l2.56|) . one may think of {—i/7i)iljx,+'4'x,- 
as the scaling limit of the Grassmann operator coupled to in the 
source term B, that is E'xj' written in terms of ip, x- From its definition, 
see Eq. (12.45 p . it is apparent that, for every finite a, such Grassmann 
operator also includes terms of the form Xx,+Xx,- or V'x.w'^f^xV'x.tj: how- 
ever, the correlations among such bilinears at distinct points vanish in 
the scaling limit. 

The truncated expectation in the last line of Eq. (l2.56|) can be graphi- 
cally represented (and explicitly computed) in terms of loop diagrams; 
correspondingly Eq.f l2.56p can be rewritten as Eq. Ol.Sp . 



2.3 The perturbed scaling limit and the temperature 
counterterm 

The last step that is convenient to perform before setting up the multiscale 
analysis that we will use to compute Eq. (l2.53p . is to properly fix the location 
of the singularity of P{d%p). Note, in fact, that the propagator associated with 
P{dil)) is singular at t = tc(0) = — 1 (and k = 0), while we know that the 
location of the singularity changes in the presence of the interaction, moving 
to t = tc{\) = tanh(/3c(A) J), where /3c(A) is the interacting inverse critical 
temperature computed in [17], which will also be derived below. Therefore, 
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it is convenient to rewrite the mass a^{0) appearing in the propagator of the 
tp field as (recalhng that tc(0) = V2 — 1) 

a JO) = ILllM = a + 47ra-V , (2.58) 
a t 

where the mass 

2 tM t - t,(A) 
" = = atJX)^~ ^ ^ ^ 

vanishes at the interacting critical point and is of order 1 with respect to a 
as a — )■ 0; we assume that a{a) ^ for every finite a > and in the scaling 
limit lima_j.o cr{a) = m*, with m* G M. On the other hand, the constant 

should be thought of as a counterterm that will be used below to fix the in- 
teracting critical temperature. The rewriting Eq. fl2.58p induces an analogous 
rewriting at k 7^ 0, i.e., a^{k) = o"(k) + A-Ka"^!/, with a{h) = a^^{cosaki + 
cos 0^2 — 2) + o". Correspondingly we decompose the inverse propagator of 

the field as C^(k) = C^(k) — 47ra~^z/cr2, where (T2 = is the second 

Pauli matrix and 

CM = Zuk)) - Q(^)C-\k)Qik) , (2.61) 

which induces the following representation for the generating function: 

S_,_(A) =kJ P,(#)P(rfx)e^(^'^'^)+^('^'^'^)+"~'"^'^"^''"^'^^ , (2.62) 

where Pa{dip) is the (normalized) gaussian integration associated to the ma- 
trix and Ma is a suitable normalization constant. Eq. fl2.62p will be the 
starting point for the multiscale analysis discussed below. But before that, 
let us discuss a few relevant symmetry properties of the fermionic action. 



2.4 Symmetries in the Grassman representation 

The integration Hicw '^^k.w'^Xk.a;, the quadratic contributions to the Grass- 
mann action -{AtiL'^Y^ Y^keVf^ i'^kCa(k.)tpy: and -(47rL^)"^ ^keD° '^-kCx(k)'?/'k, 
as well as the source and the interaction terms i3('?/', x, ^), V('?/',x, ^) and 
a~^i> J d'^il)y^(j2ipyi are each separately unchanged under any of the following 
substitutions (here we indicate Ax.j by A^, where h is identified with the 
unordered pair b = {x, x + aej)}) 
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1- ^x,<^ ioj'ip-y,,Lo, Xx,<^ -> ioJX-^,u, A A_b where, if 6 = {x,x'}, then 
—b = {—X, — x'} (parity) 

2. V^x.o; ^e^'^'^^R^-uj, X^,u> -c^e"^ixHx-<^, A withi?(xi,X2) = 
(— X2, —xi) and i?6 = -R{x, x'} = {i?x, i?x'} (diagonal reflection) 

3. i'^,uj ^V'Fx-oj, Xx,a; iXF^-uj, A Ai^'ft with F{xi,X2) = i-Xi,X2) 

and = F{x, x'} = {Fx, Fx'} (orthogonal reflection) 

4. V'x.o; ^x -t^, X^,Lo x^-u^, c ^ c*, Ah ^ Ab for all complex coefli- 
cients (complex conjugation) 

The invariance of the integration is only a matter of checking that the correct 
sign is produced, as indeed it is. 

To check the invariance of the quadratic terms, we consider the most 
general possible form which will be invariant under these transformations, as 
this will be helpful in the further analysis. A general quadratic form in the 
ip flelds is 

with C(k) = -C^(-k), that is a(k) = -a(-k), d{k) = -d{-k) and b{k) = 
— c(— k). This matrix transforms respectively as 



1. 



3. 



4. 



a(k) 6(k)\ /a(k) -c{k)\ 
c(k) d{k) J [-b(k) d(k) I ^ 



a{k) b{k)\ f-td{Rk) -c{Rk)\ 

c(k) d{k) J \ -b{Rk) ia{Rk) J ^ ''''' 

^^'^ K^)) ^ f-diFk) -c(Fk)\ ^2.66) 



c(kj d(k) J ^ V-6(Fkj -a(Fk) J 



a{k) 6(k)\ /rf*(-k) c*(-k)\ 
c(k) d(k)J \b*{~k) a*{-k)J 

Therefore, the quadratic form is invariant under these symmetries iff a(k) = 
-a(-k) = -id{Rk) = -d{Fk) = -d*{k) and b{k) = b{-k) = -c(k) = 
b{Rk) = b{Fk) = — 6*(k). Exactly the same is true of a quadratic term in 
X- We note that Co-(k) and C^^ik) as expressed in Eqs. fl2.6ip and (12.401) . as 
well as the counterterm a~^z/ J dyiip:x_a2ipx, sue indeed invariant under these 
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transformation. In passing, let us note that any quadratic form compatible 
with these symmetries is given to first order in k by 



D+(k) 



(2.68) 



with Z and a two real parameters. 

To examine the effect on V, we note that the energy density E^, (which 
is the Grassmann bilinear coupled to A^, in the source term i3(-?/;, x, A)) is 

covariant under the four symmetries above: E}^ ^ -E-fe, -Efe ^ -Em, -Eb ^ Ep^y^ 

(4) — 

E\, — > This suffices to show that the source term B is invariant under 
these transformations; furthermore we note that V is given as a polynomial 
in the E^, and A^, with real coefficients, symmetric under parity, diagonal 
reflections and orthogonal reflections, as follows from the construction in 
SecJ2.1l and therefore is also invariant. 



3 Multiscale integration 

We now want to compute H_ _(A) as expressed by Eq. fl2.62p . The strat- 
egy is to integrate the Grassmann functional integral step by step, in an 
inductive fashion. The outcome is a multiscale expansion that has been de- 
scribed in great detail in several papers, see e.g. ESI 112] for some recent 
reviews. A self-contained presentation is also described below. From now 
on C, C", c, c', . • indicate universal positive constants, whose specific values 
may change from line to line. 



3.1 The integration of the x field 

The first step simply consists in integrating out the x field, after which we 
rewrite: 

(3.1) 

where A is a shorthand for A^ and 



Recall that the label N indicates the scale of the lattice spacing, which is 
equal to a = Iq2~^ , with the unit macroscopic length (which will be set 
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equal to 1 in the following). The effective potential V^^^ on scale can be 
computed in terms of truncated expectations, namely 

|A|Fjv - log AT, + + V(^)(^, A) + A) = (3.3) 

+ E Vx (^(^' ^' ^) + ^(V^, X, A); ■ ■ ■ ; V(V^, X, A) + x, A)) 

— s times 

where £^ is defined in way analogous to Eq. fl2.57p . with P^oi^dilj) replaced 
by P{dx)- The normalization constant in the left hand side is defined 
in such a way that ^(^^(O) = V^^\0,A) = 0, while B^^\iIj,A) collects the 
contributions that are quadratic in and linear in A. Note that is a 
multilinear operator of its arguments, so that each term in the r.h.s. can be 
computed by expanding V{'ip,x,A) + B{'ip,x,A) into monomials in x and 
then by acting with on each monomial separately. The second line of 
Eq. fl3.3p can be conveniently represented graphically as in FigJBl The tree in 



N+2 

N+2 

N N+1 N N+lN+2 N N+h^ N N+VN+2 
• = • • + ^"^2 + ^ + 

Figure 6: The graphical representation of the second line of Eq.f l3.3p . 



the l.h.s., consisting of a single horizontal branch, connecting the left node 
(called the root and associated to the scale label N) with a black dot on scale 
iV + 1, represents the l.h.s. of Eq. fl3.3p . In the r.h.s., the sum of all the terms 
with s final points represents the term of order s in the r.h.s. of Eq.f l3.3p : a 
scale label N is attached to the leftmost node (the root); a scale label + 1 
is attached to the central node (corresponding to the action of where x 
is thought of as the "field on scale A^+1"); a scale label is attached to 
the s rightmost nodes ("endpoints"), which can be either "normal", in which 
case we draw them as black dots and associate them with V, or "special", in 
which case we draw them as open squares and associate them with B. We 
denote by T^^"*^ the set of trees with n normal and m special endpoints. 
Moreover, we denote the central node by Vq and the endpoints by f i, . . . , f^; 
there is a natural partial ordering on the trees from left to right, and we 
shall denote this ordering by Vo < Vi. The V and B associated with the 
endpoints consist of a sum of several different field monomials, which can be 
distinguished by assigning sets Py. of field labels /. Each field label consist of 
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a specification of tlie field type (either A, or ip, or x), and of tlie u indices of 
tlie Grassmann fields and the j indices of the A fields. We shall also denote 
by the set of field labels associated to the A and ip fields, which can be 
thought of as the fields that survive the integration P{dx)- Finally, given a 
tree r G Tj^-nm sets P^^, . . . , P^^, let P = {P^o, . . . , P^J and let Vr 

be the set spanned by P, so that the r.h.s. of Eg. (13.31) can be rewritten as 

Y: E V(^) (r, ^, A) , V(^) {r,^p,A)=J2 V^^^ (r, P) , 

n,m>0 _p7-(^) PeVr 

n+m>l ^'N;n,m 

where x„ = U^>^ U/gp^ {x(/)}, = U^>^ U/gp^ {^(/)}, and, if P^^ and P^ 
are the subsets of P^^ collecting the fields of type ip and A, respectively, 

^p.t = n ^-(/).-(/) ' = n ■ (3.5) 

Moreover, if Q,;^ = P„o I"" ^v,, noting that P^, \ Q^^ = P^^, 

4?(x.o) = ^n^f ^'Hx^J ^J(xp„V- ;XPi) , (3.6) 

where Xp^. has a definition similar to Eg. (13. 5 p and frif^^^''(x„.) is the kernel 
of the monomial labeled by P„., whose decay properties follow from Propo- 
sition [H 

The action of the truncated expectation in the r.h.s. of Eg. (13. 6 p can be 
computed in terms of a tree interpolation formula, originally due to Battle, 
Brydges and Federbush [3|, [121 [13] and re-derived in several review papers, 
see e.g. [251 126] : 

^J(xpx;-- - ;xpx) =J2^Tll9r'' f PT{t) PfG™(t) , (3.7) 

where: 

• the first sum runs over set of lines forming a spanning tree between the 
"boxes" or "clusters" vi, . . . ,Vs, i.e., T is a set of lines that becomes a 
tree if one identifies all the points in the same clusters; 

• q;t is a sign (irrelevant for the subseguent bounds); 
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• g). is a shorthand for (j/(^)(x(£) — x'(£)), where ijj{£),co'{i) and 
x(£),x'(£) are the u and x indices associated to the two ends of the 
hne i, which should be thought as being obtained from the pairing 
(contraction) of two fields Xx(^),<^(£) and x^'{e),Lo'{iy, 

• if t = {ti^i/ G [0, 1], I < < s}, then (iPr(t) is a probability measure 
with support on a set of such that tj^j' = Uj ■ Uj/ for some family of 
vectors u,; G of unit norm; 

• if 2p = X]i=i l-^j^l' then G'^~^^''^{t) is an antisymmetric {2p — 2s + 
2) X [2p — 2s + 2) matrix, whose elements are given by G^^}''^ = 

' ^'^^^^^ /' /' ^ ^eerife, fj] and f}, ff are the two field 
labels associated to the two (entering and exiting) half-lines contracted 
into i; i{f) G {l,...,s} is s.t. / G Pvi^f{, Q^^fl^) is the propagator 
associated to the line obtained by contracting the two half-lines with 
indices / and /'; 

• Pf (^^+1''^ is the Pfaffian of G^+i.^; given an antisymmetrix matrix 
Gij = —Gji, i, j = 1, . . . , 2k, its Pfaffian is defined as 

PfC = y~^(~l)'"G'^(l)7r(2) • • • G'7r(2fc-l)7r(2A:) 

7r 

= y"c?xi---c?X2fce-5S',.>^'G'i>^i , (3.8) 

where in the first line vr is a permutation of {1, . . . , 2k} and (—1)'^ is its 
parity while, in the second line, xi, . . . , X2k are Grassmanian variables. 
A well known property is that (Pf G)^ = det G. 

If s = 1 the sum over T is empty, but we can still use the Eq.f l3.7p by 
interpreting the r.h.s. as equal to if Pi is empty and equal to Pf G^+^'^(l) 
otherwise. Note that if the Pfaffian is expanded by using Eq.f l3.8p . then 
Eq. fl3.7p reduces to the usual representation of the truncated expectation 
in terms of connected Feynman diagrams. The spanning trees in Eq.f l3.7p 
guarantee the minimal connection among the clusters of fields Vi, . . . ,Vs and 
the Pfaffian can be thought of as a resummation of all the Feynman diagrams 
obtained by pairing (contracting) in all possible ways the fields x outside the 
spanning tree, with the rule that each contracted pair (xx,^, Xy,w') is replaced 
by 5'^;^/(x — y); the interpolation in t is necessary in order to avoid an over- 
counting of the diagrams. 

The reason why we prefer to use the Pfaffian expansion rather than the 
more usual expansion in connected Feynman diagrams is that the former is 
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better behaved from a combinatorial point of view: using the fact that the 
number of spanning trees in the sum Y2t r.h.s. of Eq.f l3.7p is bounded 

by s!C^^, where 2p = J2i=i l^vil (see, e.g., [251 Appendix A3. 3] for a proof of 
this fact), we find that for fixed P^^ the contribution to the kernel of ipp^i ApA 

coming from a fixed tree r G T^^^^ can be bounded as: 

E ^ / n E n i^r^^i h detc^-^^-^iiv^ (3.9) 

<\\9r''\\r' E /rfx.„nC7l^"»l|irf+^)(xJ|.||detG^+^ 

Pv-^ ^•••^Pvs i=l 



00 ' 



where || det G^^-*^'^! indicates the square root of the L°° norm of the de- 
terminant w.r.t. both the position variables x^j^ and the interpolation param- 
eter t (the square root is due to the fact that | Pf G^+^'^l = | det G^+^''^\^/^); 
moreover, in the second line, recalling that gl>^~^^'' is a shorthand for g^^^-^ (x(^)" 
x'(^)), 

\\gf^'\\i = I d^\\g\^)\\ < , (3.10) 

simply because the propagator of the x held decays exponentially on scale a = 
, as it follows from its explicit expression Eq. (12.48 p . see also Eq. (l2.5ip . 
More precisely, if 1 1 ■ 1 1 is the Hilbert-Schmidt norm, 

||'7'^(x-y)|| < C2^e-^2'^l"-yl . (3.11) 

In order to bound det G^, we use the Gram-Hadamard inequality^ stating 
that, if M is a square matrix with elements Mij of the form Mij = {Ai, Bj), 
where A^, Bj are vectors in a Hilbert space with scalar product (■, ■), then 

|detM| < JJll^ill ■ ||Ci|| . (3.12) 



where 1 1 ■ 1 1 is the norm induced by the scalar product. See [25{ Theorem A.l] 
for a proof of Eq. (13.121) . 

Let "H = M'' ® "Ho, where "Ho is the Hilbert space of the functions F : A — > 
C^ with scalar product (F, G) = E..=± / dz F*{z)G^{z), where = [F]^, 
G^ = [G]^, CO = ±, are the components of the vectors F and G. It is easy 
to verify that 

^If = U(j),i{r)gu^u),uj(f')(^(f)~^(f')) 

= (ui(/) O Bx(/),u;(/), Ui(//) (g) Cx(/'),a;(/')> , (3.13) 
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where: Uj G i = 1, . . . , s, are vectors such that tj^j/ = Uj ■ Uj/; Bx,p and 
Cx,p have components: 



[Bx,.(z)]. = J2 7 TTTi'^-.^ ' (3-14) 



so that 

^ k 0^a(k)l^+kx(k)P]'^' 

for a suitable constant C. Using the Gram-Hadamard inequahty, we find 
1 1 det G^~^'^''^\\oL^ < (const. )^'''=i 1-^^12^(3 E»=i 1^^ moreover, using Propo- 

sition [H we find that for each normal endpoint vf. 

j rfx^Jirf -*-i)(x,J| < (c|Ar^=^{i'^l^"»l}2^(2-i|P„tl-i|P^I-|P4l) . (3.16) 

Substituting these estimates into Eg. (13. 9 p we find that for fixed P^^ the 
contribution to the kernel of ipp^^Ap^^ from a given tree r G 7^^^^ can be 
estimated by 

^7i+m|j^|max{l,cn}2Af(2-||P4,|-|P4l) ^ (3-17) 

where we used the small factors |A|™'^^''^^''^'^^ to sum over the field labels, 
see ^25i Appendix A6]. Eq. (l3.17p implies the analyticity in A of the kernels 
of V(^) and of B^^\ as well as their exponential decay on scale a = . In 
particular, writing 

n>l, m>0 LO_,j i=l i=l 

(3.18) 

where the * on the sum indicates the constraint that {2n,m) ^ (2,1), we 
have 

1^ 1 rfXi ■ ■ ■ rfy™|W^£,L^^,,(Xi, . . . , X2„; yi, . . . , y^)| < C"+-2^(2-2n-„^)|^|cn 

(3.19) 

A similar expansion and similar bounds are valid for S^^\A) as well. 
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3.2 The iterative integration 

We are now left with the integration of the field. Recall that the propagator 
of the ip field is given by the inverse of Co-(k), see Eq.f l2.6ip . which has a mass 
(i.e. an inverse decay rate) proportional to cr = o"(a), which can be arbitrarily 
small. In fact, recall that a{a) 7^ Va > 0, and lim^-^o '^('3') = ^* ^ IR- In 
particular, m* = is an allowed value of the rescaled mass and it is our 
purpose to derive bounds that are uniform in m* for m* — )■ (massless 
limit). In this respect, ip is essentially a massless field or, more precisely, 
it is non-uniformly massive. Therefore, we cannot trivially integrate if) "in 
one step", as we did for x- A convenient procedure is the following. We 
define a sequence of geometrically decreasing momentum scales 2^^, with h = 
N, N—1, . . . Correspondingly we define a sequence of analytic functions fhO^) 
supported mostly around |k| ~ 2^: for instance, we can choose /Ar(k) = 1 — 
exp{-2-2(^-i)k2} and A(k) = exp{-2-2'^k2} - exp{-2-2('^-i)k2}, V/i < iV, 
so that 

1 = fM ■ (3-20) 

h<N 

The resolution of the identity Eq. (13.201) induces a rewriting of the propagator 
of ip =: ip^--'^^ as a sum of propagators concentrated on smaller and smaller 
momentum scales and an iterative procedure to compute S_ _(A). At each 
step we decompose the propagator into a sum of two propagators, the first 
approximately supported on momenta ~ 2^ (i.e. with a Fourier transform 
proportional to fh(k)), h < 0, the second approximately supported on mo- 
menta smaller than 2^, h < N. Correspondingly we rewrite the Grassmann 
field as a sum of two independent fields: = ip^^') + -ipi^^-^) and we inte- 

grate out the field ijj^'^^ in the same way as we did for x- The result is that, 
for any h < N, we can rewrite Eq. (l3.1l) as 

S_ _(A) = e^MFh+SW(A) j p^^ ^^^^ (^^^(</i)^)gV('^)(^(<'^),A)+BW(^(<'>),A) ^ 

(3.21) 

where Fh, Zh,ah,V^^\B^^'i will be defined recursively, X/i(k) = Yjk<hfk{^) 
and Pxh,Zh,o-hi'^''P^~'^^) gaussian integration with propagator (recall the 

definition D"^ = a~^{i sin aki ± sin 0/^2)) 

„(<.)^^) _ 2vr^ e-'--X.(k) /Z^+(k) ^a, \ 
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li h = N, then: Z^r = 1, ctjv = cr(ct) = cr and 

(3.23) 

In the following steps, the effective potential V^^^ and the effective source term 
B^'^^ will be shown to have the following structure: B^'^\ip,A) is quadratic 
in ip and linear in A, while V*-^^ admits an expansion analogous to Eq. fl3.18p 



* „ 2n m 



?i,m>0 oJ_,j i=l i=l 

?i+m>l ~ 

(3.24) 

where, as we will see below, 

1 rfxi ■ ■ ■ rfy„|iyi„\^,^(xi, . . . , X2„; yi, . . . , y„)| < C"+-2'^(2-2n-H|Ap 

(3.25) 

The iteration continues until the scale h = h„ := [logj (^(a)] is reached. At 
that point, the left-over propagator, g^-^'^^ is massive on the "right scale" (i.e. 
on the very same scale 2'^'^), so that the associated degrees of freedom can 
be integrated in one step. The result is the desired generating function, from 
which we can finally compute the multi-point energy correlation functions. 

3.3 Localization and renormalization 

In order to inductively prove Eq. fl3.21l) we write 

V(^) {ijj, A) = rv^'^) i^) + 7^V('^) i^, A) , (3.26) 

where, if we think of the kernel Wl'l! , s as a 2 x 2 matrix with matrix 
indices ui,U2, 

CV^^\^) = (3.27) 



L2 



+ Yl / ^^^1 ■ ■ ■ t^X4W^4,0;^(xi, X2, X3, X4)^gt\ ^gi^.^g^^^gil , 



UJi ,UJ2 ,i^4 



Here Vq (resp. Vi) is an operator that extracts the order (resp. order 1) in 
Wk}k>h from the kernel which it acts on, and da(k) := (sin a/ci, sin 0/^2) • 

Remarks. 
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1. L will be called the localization operator, which should be thought of 
as the linear operator extracting from the effective potential V*^'*^ its 
local (singular) part, while TZ will be called the renormalization opera- 
tor, which is the linear operator extracting from V'^'^'' its regular part. 
Note that the action of TZ on the kernels quadratic in ip is equiva- 
lent (via the use of the remainder's formula in the Taylor's expansion) 
to the action of a second order differential operator of the type k^5^ 
or Xlfc>/i'^fc^CTfc (^^^ ^ suitable combination of the two). Dimensionally, 
the differential operator k^(9^ (or, cquivalently, its x-space counterpart, 
which is of the form x^S^) behaves as 2'^^'^^~'^^\ where hi is the scale 
of k and /i2 the scale of the propagator that acts on; moreover, 
the iterative integration is set up in such a way that by construction 
hi < h2, so that 2^*^''i^^2) jg ^ dimensional gain, sufficient to regularize 
the quadratic kernels of the effective potential. Similarly, the action of 
the operator Ylk>h '^l^lk '^'^ kernels quadratic in ip is dimensionally 
equivalent to a multiplication by cr^2^^'' ~ (^.£k.Y2'^^^''~^^; as we will 

see, Gh — (T — a{a) at all scales, so that the action of J2k>h'^k'^lk 
dimensionally equivalent to a multiphcation by the gain factor 2'^'^^""^^ 
which is enough to regularize the quadratic kernels of the effective po- 
tential. A similar discussion is valid for the action of TZ on the quartic 
kernels. 

2. A key fact which makes the theory at hand treatable (and asymptot- 
ically free) is that the quartic term in the second line is zero "by the 
Pauli principle", i.e., simply by the Grassmann rule V^^.a; — 0- ^Si.ct, 
note that at least two of the four ui indices must be equal among each 
other. Therefore the integrand in the second line is identically zero. 
This property can be diagramatically interpreted by saying that the 
fermionic nature of the theory automatically renormalizes the four-field 
interaction, which is dimensionally marginal (see below) but effectively 
irrelevant thanks to the cancellation that we just mentioned. 

Similarly, we decompose the source term as 

S^*) (V', A) = CB^^^ {i^,A)+ UB^^^ {^,A) , (3.28) 

with 

£^W(^,A) = -i (4,i + ^,2)^-k-pWff(0)4, (3.29) 

k,pex>M 

where = J dxe^P^Axj. As we will see below, CB'^^'^ is a marginal operator 
in the Renormalization Group sense. 
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The symmetries of the theory, which are described in Sec J2.4l and are 
preserved by the iterative integration procedure, imply that the kernel in the 
second line of Eq. fl3.27p has a structure analogous to Eq. fl2.68p ; we also choose 
to separate the mass term (i.e. the constant in the off-diagonal elements of the 
matrix) in two parts, one proportional to the bare mass and one independent 
of it, namely: 



k ^ ' 



=: £oV(")(^) + 2VF,(7A), (3.30) 

where F^{ip) := f d'^ipy^a2'ipyi and Ch,Sh,i^h are suitable real constants, such 
that Sh is linear in {ak}k>h and Uh is independent of {<7k}k>h- Note that the 
counterterm on scale h = N is defined so that 

z/jv := z/ + 2-^-'TT[cr2VoWiTW] ■ (3-31) 

We stress once again that the local part of V*-^-* is purely quadratic in tp: the 
quartic term, a priori present in CV^'^^ is zero thanks to the "Pauli principle" , 
in the sense of the Remark 2 above. Similarly, the local part of the source 
term can be written as 

jr^W ^ / + A,,2)V'xa2^x . (3.32) 

2vr J 

for a suitable real constant z\^\ with Z^^^ = 1. 

Once that the above definitions are given, we can describe our iterative 
integration procedure for h < N. We start from Eq. fl3.2ip . which is induc- 
tively assumed to be valid at the h-th step, and we prove the validity of the 
representation for h — 1. We rewrite Eq. (13.21 p as 



(3.33) 



Next we include CqV^'^^ in the fermionic integration, so obtaining 

.g7^v('•)(v<-''',^)+B<'''(V'*-''^^) ^3 

where 

Zh-i{k) = Zh + ChXh(M) , Zh-i{k)ah-i{k) = ZhOh + ShXh{k) , (3.35) 
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and Ch is a constant fixed in such a way that / Px^,z^_„TT,,_,{d^^-^y) = 1. Now 
we can perform the integration of the tjj^'^^ field. We rewrite the Grassmann 
field as a sum of two independent Grassmann fields ijj^-'^ + tjj^'^^ and 

correspondingly, if we let Z^-i := Zh-i{0) and cr^-i := Wh-i{0), we rewrite 
Eq.([S31 as 



where 



A(k) 



Xh{k) 




Zh-i{k) 





Xh-i(k) 



-1 • 



Zh-i 

A(k)?,_i(k) = [A(k) + x/.-i(k)]a,_i(k) - x/.-i(k)a. 
The single scale propagator is 

^ ^ Zn-,L\^jDa{k)\^ + al_,{k)\-ian-M ^-(k) 

(3.36) 

A key remark is that, if \Zh — 1| < (const.) |A| and |cr/j — cr| < (const.) |o"A| 
(as we shall inductively prove below) and for all scales larger than h^, the 
propagator g^^^ satisfies a bound analogous to Eg. (13. 11 1) : 

\\g^''\^-y)\ \ < C2V'=''l"-yl . (3.37) 

Therefore, we can integrate the field on scale h in the same way as we did for 
X, after which we are left with an integral involving the fields and 
the new effective interaction, defined as 

= I Pj , ~ (rf^W)e^<'^'(^'^''-^'+^<'^^^) , 

where 

V^^\^) = 2^VhFM + CB^^\i^. A) + 7^V(^)(^> A) + TZB^''\^, A) , 

and V'-^'^\0,A) = 0, so that = Fh + eh + Chand ^(^-^^(A) =5W(A) + 
Sh-i{A), while B^'^~'^\ijj , A) collects all the terms quadratic in tjj and linear 
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in A. The potential V^'^ is of the form Eq. (l3.24p as one sees by using the 
identity 

\A\eh + SH^iiA) + V^''-'\i^^^''-'\A) + B^''-'\i^^^''-'\A) = (3.39) 

= y l-g^(V^'^\^p(^'^-') + A);...; V^'^H^^"""'^ + A)) , 
^ ' 

— n times 

where £^ is the truncated expectation w.r.t. the propagator g^'^\ 

The iteration continues up to scale h„. At that point, we integrate in one 
step all the remaining degrees of freedom, by taking advantage of the fact 
that 

- y)|| < C'2'^<^e-'=2'"l"-yl , (3.40) 

where by definition 2^"^ is of the same order as a = a (a) ^ 0. After 
the integration of ip^-^'^^ we are left with the desired generating function, 
logS_,_(A). 

Note that the above procedure allows us to write the effective constants 
{Zfi, ah, i^hi Z^^) with h < N, in terms of {Z^, cr^, z/^, Z^"*) with h < k < N: 

Zh-, = Zh + (3^, ^ = l + /3-, (3.41) 
z/,„i = 2z/, + /3^ ^ = i + /3fi, (3.42) 

where /3f = /3f{{Zh, at, J^h, zj^^), • • • , {Zn, o^n, vn, Z^^^)) is the so-called Beta 
function. 

The effective constants, sometimes called the running coupling constants 
measure the strength of the local part of the effective potential. The key 
point, to be proved below, is that, if we assume that the running coupling 
constants stay close to their bare values at scale A^, then the expansion for the 
effective potential induced by the iterative procedure above, is well defined 
and analytic in the sequence of running coupling constants. Moreover, under 
the same assumptions, the beta function itself is well defined and analytic 
in the running coupling constants: this allows us to study the evolution of 
these effective constants under the dynamical system defined by We will 
see that, uniformly in 

Z, = 1+0(A), Z« = 1+0(A), Vh = 0{X), (t, = (t(1+0(A)), 

(3.43) 

where a = a (a). The boundedness of the flow of the running coupling con- 
stants, Eq. fl3.43p . is one of the key ingredients that allow us to prove the 
analyticity of the theory in the bare couphng constants, see below. 
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3.4 Tree expansion for the effective potential 

In order to prove the analyticity properties of the effective potential an- 
nounced in the previous section, we first need to prove that the kernels of 
the effective potential, as obtained from the iterative procedure described 
above, are expressed by absolutely convergent series in the sequence of the 
running coupling constants, provided that these are assumed to be close to 
their bare values at scale h = N. Then we will show that this assump- 
tion is justified, by proving that their flow under the beta function remains 
bounded. In any case, as a first step towards the full control of the theory, 
we need a more explicit representation of the kernels of the effective potential 
and a systematic and efficient way to bound them. The idea is to systemat- 
ically represent the action of S"^ in Eq. fl3.39p in a way analogous to the one 
described in Section 13.11 for the integration of the x field. By iterating the 
graphical equation Figj6] we obtain a tree expansion for the effective poten- 
tial, which has been first introduced by G. Gallavotti and F. Nicolo in [23] 
and since then it has been described in detail in several papers that make use 
of constructive renormalization group methods, see e.g. [21j and the more re- 
cent reviews [251 [261 H2] . The main features of this expansion (to be referred 
to as the expansion in Gallavotti-Nicolo trees, or in GN trees for short) are 
described below. 



^^0 




h h+1 N+2 

Figure 7: A tree r G Tj^^nm with n = 7 and m = 2: the root is on scale h 
and the endpoints are on scales < + 2. 



1. Let us consider the family of all trees which can be constructed by 
joining a point r, the root, with an ordered set of n + m > 1 points. 
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the endpoints of the unlabeled tree, so that r is not a branching point. 
The endpoints can be of two types, either normal or special, the first 
drawn as black dots, the second as open squares, see Fig|71 n and m 
indicate the number of normal and special endpoints, respectively. The 
branching points will be called the non trivial vertices. The unlabeled 
trees are partially ordered from the root to the endpoints in the natural 
way; we shall use the symbol < to denote the partial order. Two 
unlabeled trees are identified if they can be superposed by a suitable 
continuous deformation, so that the endpoints with the same index 
coincide. We shall also consider the labelled trees (to be called simply 
trees in the following); they are defined by associating some labels with 
the unlabelled trees, as explained in the following items. 

2. We associate a label < h < N with the root and we denote by T^^l ^ 
the corresponding set of labeled trees with n normal and m special 
endpoints. Moreover, we introduce a family of vertical lines, labeled 
by an integer taking values in [h,N + 2], and we represent any tree 
r G TJ^'^^m that, if v is an endpoint, it is contained in a vertical hne 
with index h + l<hv<N + 2, while if it is a non trivial vertex, it is 
contained in a vertical line with index h < hy < N + 1, to he called the 
scale of f ; the root r is on the line with index h. In general, the tree will 
intersect the vertical lines in set of points different from the root, the 
endpoints and the branching points; these points will be called trivial 
vertices. The set of the vertices will be the union of the endpoints, of 
the trivial vertices and of the non trivial vertices; note that the root 
is not a vertex. Every vertex f of a tree will be associated to its scale 
label h^, defined, as above, as the label of the vertical line whom v 
belongs to. Note that, if Vi and V2 are two vertices and Vi < V2, then 
hu-^ <^ • 

3. There is only one vertex immediately following the root, called Vq and 
with scale label equal to h + 1. 

4. Given a vertex of r G T^]^^ that is not an endpoint, we can con- 
sider the subtrees of r with root f , which correspond to the connected 
components of the restriction of r to the vertices w > v. If a subtree 
with root V contains only v and one endpoint on scale h^ + l, it will be 
called a trivial subtree. 

5. With each normal (resp. special) endpoint v on scale A^+2 we associate 
a factor V{^^^^\ x, A) (resp. B{^p^^^\x, A)); here x =■ i^^^^^^ should 
be thought of as the field on scale + 1. With the endpoints on 
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scale < N + 1 we associate a factor 2^''~'''Vh^_iFy{il)^-^''~^^) if the 
endpoint is normal or a factor CB^'^^'^^ if the endpoint is special. The 
vertex v' immediately preceding an endpoint v on scale < N + 1 is 
necessarily non trivial. Note that none of the endpoint on scale < + 1 
is associated with a quartic operator in the ip fields: this is due to the 
cancellation mentioned in Remark 2 at the beginning of Sec J3.3[ 

In terms of these trees, the effective potential V^^-* can be written as 



where, if Vq is the first vertex of r, if ri, . . . , (s = Sy^) are the subtrees 
of r with root vq, and if S]^_^_i is the truncated expectation associated to the 
propagator ^f'^'^^. 



(3.45) 

and F^''+'\ri,V^(^'^+i),A): 



• is equal to TZV^'^'^^\Ti,'ip^-'^^^\ A)) (where TZ is the linear operator 
induced by the definitions Eqs. ( 13^ . (13:271) .( 13:281) .( 13:291) ) if r, is non 



• is equal to 2^+^z//i+iFy(^/'^-^"'"'^^) if Tj is trivial, h < N and the endpoint 
of Ti is normal; 

• is equal to CB^'^'^^\'iJj^-'^'^^\ A) if Tj is trivial, h < N and the endpoint 
of Tj is normal; 



• is equal to V(V^(-^+^\ A) (resp. i3(^(^^+i\ A)) if r, is trivial, h = N 



and the endpoint of Xj is normal (resp. special). 

Note that, by the definition of V^^\ it is apparent that all the nodes that are 
not endpoints are associated with the action of the renormalization operator 
TZ. The local (singular) parts only appear in the contributions associated with 
the endpoints: in this sense, the singular parts "do not accumulate" . Using 
its inductive definition Eqs. (l3.44p - (l3.45p and the Pfaffian representation for 
the truncated expectations, the right hand side of Eq. ( 13.44]) can be put in 
a form similar to the one derived in Section 13.11 for the integration of the 
X field. To do that, we further expand V^'^\t,^P^-'^\ A) by distinguishing 




(3.44) 



n,m>0 pT-Ci) 

n+m>l ^'jV;n,m 



trivial; 
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the different contributions arising from the choices of the monomials in the 
factors V and B associated with the endpoints on scale + 2, as well as the 
scale at which each field in these monomials is contracted. To this purpose, 
we introduce a few more definitions, which generalize the definition of P„ in 
Sec J3.1l 

We introduce a field label f to distinguish the field variables appearing in 
the monomials associated with the endpoints; the set of field labels associated 
with the endpoint v will be called /„; if f is not an endpoint, we shall call 

the set of field labels associated with the endpoints following the vertex 
V. Note that every field can be either of type A or ip: correspondingly, we 
denote by and the set of field labels of type A and ip, respectively, 
associated with v. Furthermore, we denote by x(/) the space-time point of 
the field variable with label /; if / G I^, we denote by the j index of 
the external field with label /; if / G I^, we denote by the u index of 
the Grassmann field with label /. 

We associate with any vertex v of the tree a subset of 1^, the external 
fields of v; we further denote by and P^ the subsets of Py of fields of 
type A and V, respectively (of course, P^ r\ Pf = % and P^ U Pf = P^). 
These subsets must satisfy various constraints. First of all, if v is not an 
endpoint and vi,...,Vs^ are the s^, > 1 vertices immediately following it, 
then Py C UjP^j if v is an endpoint, Py = ly. If f is not an endpoint, we 
shall denote by Q^^ the intersection of P^ and P^. ; this definition implies that 
Py = UiQy.. The union of the subsets \ Q^^ is, by definition, the set of 
the internal fields of v, and is non empty if s„ > 1. Similar definitions are 
valid for Q^, Q^, etc. Note that P^ = for all v, simply because A is, 
by definition, an external field. Given r G T^'^lm' there are many possible 
choices of the subsets P^,, f G r, compatible with all the constraints. We 
shall denote by Vr the family of all these choices and by P the elements of 

Let us note that the resulting expansion for the effective potential, as 
compared to the expansion for V^^^ described in Sec. 13.11 has an important 
difference related to the iterative action of the TZ operator on the nodes of 
r that are not endpoints; as observed in the first Remark at the beginning 
of Sec J3.3l this is equivalent to the action of a suitable differential operator, 
whose precise form has been discussed in several papers on the subject, see 
e.g. [SmnillT]; fortunately, in our context we do not need to describe the form 
of the interpolation operator exactly, but only some of its general properties, 
discussed in the following. 
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3.4.1 The non-renormalized expansion 

Let us start with describing the basic dimensional bounds of the effective 
potential, temporarily neglecting the action of the renormalization operator. 
More precisely, let us temporarily pretend that the action of TZ on the nodes 
of T that are not endpoints, induced by definition Eq. fl3.45p and by the first 
item immediately following it, is replaced by the identity. Then the result of 
the iteration would lead to the following relation: 

Vf)(r,V^(^'^), A) =J2J2f ^^-oW^;,P,T(x.o)<f , (3.46) 

where Xt,^ is the set of integration variables associated with r and T = 
Ui. not ep '^v is the union of the spanning trees associated with all the nodes 
that are not endpoints in r. Moreover, W*p rj, is given by 

W^;,p,t(x.o) = (3.47) 

= [ n { n ^ / ^^^"(t^) p^G'-^'i^^) [ n ^j'^^] } ' 

V e.p. V not e.p. ^ i&Tv 

which can be thought of as the multiscale version of the expansion derived 
in Section 13.11 Here s„ indicates the number of nodes immediately following 
V on r. 

Eq. fl3.47p can be bounded in a way analogous to Eq.f l3.9p : roughly speak- 
ing the rationale is that each gf^"^ is replaced by its Li norm, which dimen- 
sionally is proportional to 2"^", and each Pf G''^"'"^"(tt,) by the square root of 
its Loo norm, which dimensionally is proportional to 2 ""-a ^i=i l^^il-si-^i' l-^v+i-) _ 
The result is 



i- J d^jW:^^^A^J\ < [ n (C|A|)— {i'^l^"l>] • (3.48 

P ti normal e.p. 



(2-i|P„'*|-|P„^|)' 



„ 2 1-"" 

i;not e.p. " v e.p, 

where in the last factor v' indicates the node immediately preceding v on 
r. Note that the bound is valid provided the running coupling constants 
are close to their bare values. The big product in the second line can be 
conveniently reorganized by using the identities 



h^{sy - 1) = h{n + m - 1) + ^ {h^ - h^/) {n{v) + m{v) - 1) , 

not e.p. tinot e.p. 



V not e.p. i=l v not e.p. 
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where n{v) (resp. m(f)) is the number of normal (resp. special) endpoints 
following V on r. In the second line, |/^| = \P^\, so that, plugging 

these identities into Eq.f l3.48p and using also the fact that J2v ep 1-^? 



A\ 



V e.p. I V 



l^^^l, gives 



normal e.p. 

.2M2-||p4,|-|p4|) J j-j- 2('^"-'*-')(2-|l^.*l-l^'„^l)" . (3.49) 

Dnot e.p. 

We call djj = —2 + — \P^\ the scaling dimension of f , depending on the 
number of the external fields of v. If < for any v one can sum over r, P, T 
obtaining convergence for A small enough; however dy <0 when (|-P^|, \P^\) 
is equal either to (2,0), in which case the scaling dimension is 1 (relevant 
operator), or to (4,0) or to (2, 1), in which cases the scaling dimension is 
(marginal operator). The apparent divergences associated with the presence 
of nodes with these special fields configurations is cured by the action of the 
regularization operator TZ, as discussed in the following. 



3.4.2 The renormalized expansion 

Of course, the expansion for V*^^-* described above does not lead to Eq. (13.461) . 
but rather to a different expansion, which is similar in many respects to 
Eq. (13.461) . modulo the presence of a certain number of interpolation operators 
arising from the action of TZ on the nodes of r that are not endpoints. The 
precise definition of these interpolation operators is rather technical, see e.g. 
Section 3 of ^Jj. The outcome can be expressed as 

nV^^\T,^^^''\A) = (3.50) 



E E E / ^x.oiy.,p,T,,(x.j[ n 

"o^'n T-'r-Ti a^TD J _ _ 



5(/)./,(<^) 

''x^(/),a;(/) 



A 



pA 

"0 



where: Bt is a set of indices which allows to distinguish the different terms 
produced by the non trivial TZ operations; x^(/) is a coordinate obtained 
by interpolating two points in x^,;,, in a suitable way depending on /3; g^(/) 
is a nonnegative integer < 2; jpi^f) = 1,2 and (9J is a suitable differential 
operator, dimensionally equivalent to 9J (see [TT] for a precise definition); 
Wr,p,T,/3 is given by (defining Si = \ — and 1S2 = 1 — Po — "^i, as in [27] : 
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recall that Vq and Vi were defined at the beginning of Sec J3.3p : 



V e.p. 

} 



V not e.p. 



(3.51) 



where: (ij^(£)(x^ — x^) is a suitable interpolation operators, dimensionally 
equivalent to |x£ — x^| (and similarly (ij^(^)(x„j3)), is dimensionally equiva- 
lent to the distance operator among two distinct points in x^); the indices 
bpiv), hp{l), q^Ul) and qpUl) are nonnegative integers < 2; ip{v), jpUl)^ 
j/sifi) and j/3(i') can be either 1 or 2; ^/3(f) and ii3{i) can be either 1 or 
2; //3(f) and //?(/) can be either or 1; C/^i^v), c^(f), C'/3(Z) and C/3(/) can 
be either 0,1 and max{C/3(f) + c^(f),C/3(/) + c^(0) ^ Ij G*^"' "(^f) is ob- 
tained from G'^'"'^'' (ty) by substituting the element ti(^f-j ii^f/^g^^^^f, f) with 

'<(fl.-</')a:;;,'9i',%'9""'(/,/'). 

It would be very difficult to give a precise description of the various con- 
tributions of the sum over Bt'- however, we only need to know that the 
number of "zeros" dj^(^i){'Xi — x^), the number of "derivatives" dj^(^f) and the 
scale of the propagators they act on are, by construction, such that all the 
potentially dangerous nodes v that are not e.p. (i.e. the nodes v that are 
not e.p. and with {Pf,P^) equal either to (2,0), or (4,0), or to (2, 1)) gain 
a "scale jump" 2"(^''"^"')^(-^''), where 



n if (inl,in^i) = (4,o), 

2 if (|P.1,|n^|) = (2,0), 

1 if (|n1,|n^|) = (2,l), 

otherwise 



(3.52) 



which is enough to cure the apparent divergences present in Eq. fl3.49p . For 
more details, see e.g. [HI |25l [27] , see also [281 Sec. 2. 2] for a recent pedagogical 
description of this point. 

In conclusion, the kernels of the effective potential are bounded in a way 
completely analogous to Eq. (13.491) . modulo the improved dimensional factors 
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induced by the action of the TZ operators: 

\max{l,c|P^|} 



/3G-Bt 11 normal e.p. 

.2M2-|I<|-|p4I) J Y[ 2(''"-'*''')(2~5l-P-'*l-l^.''l-<-P-))" , (3.53) 



Dnot e.p. 



which is vahd provided that the running coupling constants stay close to their 
values on scale h = N, for all scales between h and A^. Under this assump- 
tion, Eq. fl3.53p implies the analyticity of the kernels of V*-^^ and B^^^ and the 
decay bounds Eg. (13.251) . 

An immediate corollary of the bound Eq. (13.531) is that contributions from 
trees r G TJ^^^ ^ with a vertex v on scale hy = k > h admit an improved 
bound with respect to Eg. (13.251) . with an extra dimensional factor 2^^'^~^\ 
< 9 < 1, which can be thought of as a dimensional gain with respect to 
the "basic" dimensional bound in Eq. (13.251) . This improved bound is usually 
referred to as the short memory property (i.e., long trees are exponentially 
suppressed); it is due to the fact that the renormalized scaling dimensions 
djj = 2—^\P^\ — \P^\—z{Py) in Eq. (l3.53p are all negative, and can be obtained 
by taking a fraction of the factors 2^'^""'*"')'^" associated to the branches of 
the tree r on the path connecting the vertex on scale k to the one on scale 
h. 

Under the same assumptions, the beta function itself, is analytic 
and dimensionally bounded by a constant independent of h. Moreover, the 
contributions to it from trees that have at least one node on scale k > h 
is dimensional bounded proportionally to 2^^^''~^\ with < < 1. It is 
remarkable that thanks to these bounds, the dynamical system induced by 
the beta function can be fully studied and shown to lead to a bounded flow 
of the running coupling constants, as discussed in the next section. 



3.5 The flow of the running coupUng constants 

As announced above, the flow of the running coupling constants is controlled 
by the equations: 

z/,_i = 2z/^ + /3^ ^ = i + /3fi, (3.54) 

Zh 

where /3f = Pf{{Zh, ah, i^h, ^I^^), • • • , {Zn, ctn, ^n, ^^^)) is an analytic func- 
tion of its argument, with an analyticity domain bounded by: \Zk — 1| + 
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— 1| + |z/fc| + 1^ — l| < eo, for all /i < k < N, where £0 is a suitable 
(small) positive constant. 

Let us start by studying the flow of Zh and u^. By construction, their beta 
functions, and /3^, are independent of zj^'' and ah and, in particular, they 
are independent of the infrared cutoff scale h„, for all h > hm- Therefore, the 
two coupled equations for Zh and Uh can be naturally studied for all h < N, 
without any infrared cutoff. Note also that both (3^ and (3^ can be expressed 
as sums over GN trees with at least two endpoints, and with at least one 
endpoint on scale + 2, the reason being that the local part of the trees 
with only normal endpoints of scale < + 1 is zero by the support properties 
of the single-scale propagators that enter the definition of /3j^. Therefore, by 
the short memory property, \/3h\,\(^h\ — Ce|A|2^'^^^^^ 

In order to solve the two coupled equations for Zh and Uh we use a fixed 
point argument, following the same strategy as [271 Appendix A5]. Let 
^N;K,e be the space of sequences u = {i'h}h<N such that < K\\\2~^^'^~-^\ 
\fh < N] we shall think ^fftN;K,e as a Banach space with norm || ■ H^, where 
llz^lle = sup^<jv 2^^^~^''- Given u G ^n;K£, we first solve the equation for 
Zh, Zh-i = Zh + f3h , which leads to a bounded and exponentially convergent 
flow. In fact, using the fact that \(3^\ < C\\\2^^^~^\ it is straightforward to 
check inductively that, if Z^ = 1, 

\ZhiK) - 1| < C|A| , \Zh{K) - Z^Uk)\ < 01X12"^ (3.55) 

uniformly in uforuE ^N]K,9- Once that Zh = Zh{v) is fixed as a function of 
the beta function for Vh can be thought of as function of u_ only, and we shall 
write Phi^)- The goal is to fix un in a smart way, so that the flow generated 
by /5^(z^) with initial datum generates a sequence in DyiN-K,e- We note that 
the solution to the flow equation induced by the beta function with initial 
condition lim/i_j._oo zz/^ = is a fixed point of the map T : Dyij\f.K,e ^N;K,e 
defined by 

(Tz£). = -5^2^--^-i/3J(z.). (3.56) 

The fact that, for K sufficiently large, T is a map from dyt]s[-K,e to itself is 
a simple consequence of the bound < C|A|2^(^-^). More interestingly, 
we can prove that T is a contraction map on dytN;K,e- In fact, using the 
expansion in GN tress, the short memory property and the fact that the 
trees contributing to (3^ have at least one endpoint on scale + 2, we find 
that if z/, I/' G TlN-K,e, then 

mi^) - (3hiK)\ < C|A|5^|z/,-z/^|2^('^-^) < C"|A|2^('^-^)||z/-z/'||, , (3.57) 

k>h 
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which imphes, if plugged back into Eg. (13.56 p . that 



\\Tu-Ti/\\g < c"|A|||z/-z/'||(, , (3.58) 

i.e. T is a contraction for |A| sufficiently small. Then the Banach fixed point 
theorem shows that T has a unique fixed point u* in ^fft^-Kfi, such that 

^N = i^*N = -Yl 2'"'^"'/3;(i^*) , (3.59) 

j<N 

which is of order A, as it should be. In fact, plugging Eq. (l3.59p into Eq. (l3.30p 
and using the relation of u with the perturbed critical temperature, Eq. (l2.60p . 
shows that the distance between the interacting and the free critical temper- 
atures is of order A, as expected. A close examination of the proof also 
shows that the resulting u is exactly independent of a (rather than being 
"just" bounded by C|A|, independently of a). Similarly, the limiting value of 
Zh = Zhijj^) as h ^ — oo, to be called simply Z_oo(A), is exactly independent 
of a. 

Once that the fiows of and have been controlled, the fiows of ah 
and zi^^ can be studied easily, in the same fashion as the fiow of Z^ at 

Z 1 

jy G ^N;K,e fixed. First of all, note that by construction both and 
are independent of ah and, in particular, they are independent of the infrared 
cutoff scale h^j, for all h > h„. Therefore, also the two equations for 

^(1) 

and — ^ can be naturally studied for all h < N, without any infrared cutoff. 

Moreover, f3'^, in complete analogy with the beta functions for Zh and Uh, can 
be expressed as sums over GN trees with at least two endpoints, and with 

Z 1 

at least one endpoint on scale N + 2; similarly, is a sum over GN trees 
with exactly one special endpoint (on an arbitrary scale) and at least one 
normal endpoint on scale + 2. Therefore, by the short memory property, 
I, \Ph^\ < Ce\X\2^^''-^\ from which we get the analog of Eq. fl33HD (recall 
that (Tat = 0" and Z^^ = 1): 

\ZIP - 1| < C|A| , - Zi'il < C|A|2^'^ , (3.60) 

I— -1|<C|A|, I— -Z*(A)| <C|A|2^\ (3.61) 
' cr ' 'a 

where Z^l^{X) and Z*{X) are two analytic functions of A, analytically close 
to 1. As already observed for Z__oo(A), these two functions are independent 
of a. 
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4 The renormalized correlation functions 



We are left with computing the renormahzed m-points energy correlation 
functions and proving the main results stated in Section [1] and summarized 
in Theorems 11.11 and 11.21 We closely follow |71 Section 2.3]. Recall that 
the m-point correlation in the thermodynamic limit and at distinct points 
Xi, . . . ,Xm, can be obtained by performing the functional derivative of the 
generating function H_^_(A) constructed in the previous sections (here /3 = 
/3(a), fixed so that 2(t — tc{0))/t = aa + 47rz/, where t = tanh(/3(a) J), u is 
the temperature counterterm fixed as explained in the previous section and 
cr = a{a) is an a priori prescribed mass, such that lim^^o "^(c^) = f^* ^ ^)'- 



=0 

(4.1) 

Performing this functional derivative is just a way of extracting from the 
generating function the kernel associated with the monomial ^xi.ji ■ ■ ■ ^x^jv™ 
in logS_„(A) = J2h=h -i<Sh{A), where Sh{A) with h < N is defined by 
Eq. flQQj) and Sn{A) := S^^^A), see Eq.([33]). The function Sh{A) can be 
written in a way similar to Eq. fl3.24p : 



/m 
dxi ■ ■ ■ cix^5i,^.^^(xi, . . . , x^) JJ A^^j. , 

•rri^i- Jl, ■■■,3m i=l 



(4.2) 



so that 

N 

(exi,ii; ■ ■ ■ ; e^mjmfp =m\ Si^^li^i, . . . , x„) (4.3) 

h=hcr — l 

and the thermodynamic limit L — )• oo is understood. Note the combinatorial 
factor ml arising from the action of the derivatives w.r.t. A^.j- on Sh{A). 

As explained in the previous sections, S'^?^ (xi, . . . ,Xm) can be expressed by 

a sum over trees r G T^^l^l with niQ < m special endpoints and n > 
normal end-points, and over a suitable set of field labels P & Vr', ttiq < m 
corresponds to the case where there is at least one normal endpoint v on 
scale + 2 such that ^ 0. In formulae, we can write: 

(S^xiji! ■ ■ ■ i £^Xm,jm)/3 ~ (4-4) 

N m 



) Jl) ■ ■ ■ ) ^m; 3 

h=h^~l n.>0 ■mo=0 PePr: 



m j } 



58 



h K N + 2 

Figure 8: Example of a tree r G TJ^'*^ appearing in the expansion for the 
m points correlation function, with n = 3 and mo = m = 3. The subtree 
r* associated with r and with a choice of P such that \Pvo\ = l-P^I = 3 is 
highlighted. 



where 5't-^p(xi, ji; • • • ; x^, jm) can be represented as in Eqs. fl3.50p - fl3.5ip . and 
bounded in a way similar to Eq. fl333|) . Given r G T^^i,™^ and P e Vr, let 
£'*p be the set of endpoints v in r such that ^ (i.e. E*p contains all 
the special endpoints of r plus the normal endpoints v on scale N + 2 such 
that \P^\ > 1); let us denote by r* the minimal subtree of r connecting all 
the endpoints in E*p. For each f G r, if := \P^\, let s* the number of 
vertices immediately following v with m„ > 1. Moreover, let Vnt{T*) be the 
set of vertices in r* with s* > 1, which are the branching points of r*. For 
future reference, we also define Vq to be the leftmost vertex on r* and its 
scale. See FigjHl 

Given these definitions and recalling that dy(P^) = 2 — — \P^\ — 

z{Pj;), we can write the bound for 5't-^p(xi, ji; ■ ■ ■ ; x.^, jm) as 

|^.,p(xi,ji;-- - ;x^,j^)| <C"+™2^(2-")[ JJ 22(<-i)'^V-'=2''"^^] ■ (4.5) 

( |A|'^I-^"'I, if V is normal, = N + 2 
■[ n 2(''--''"')'^"(^")] J] ^ |A|2^('^-^), if is normal, /i, < iV + 2 

''note.p. ?; e.p. [^1^ if ^; ig special 

This estimate is very similar to the bound Eq. (13.531) for the renormalized 
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kernels of the effective potential, the most important difference consisting 
in the product ni;ev'„t(r*) 2^('*^~^)''"e~''^'''''^", where 5^ = 5i,(xi, . . . ,Xm) is the 
tree distance of the set x* := UjgpA{x(/)}, i.e. the length of the shortest 
tree graph on connecting the points of x*. In this product, the factors 
22(«S-i)'»« take into account the dimensional gain coming from the fact that 
we are not integrating over the space labels x, of the external fields (the gain 
is meant in comparison with Eg. (13.531) where, on the contrary, we integrated 
over all the field variables); moreover, the factors e~'^^ come from the 
decaying factors e~^'^ ^|x(£)-x'{£)| associated with the propagators g'^'" on the 
portion of spanning trees inside the cluster v, that is from the lines i E T^j, 
w > V, see [3 Section 2.3] for a few more details about how to extract these 
factors starting from the bounds on the kernels of the effective potential. 
Moreover, note that in the second line we estimated the contribution from 
each endpoint of type u by (const.) |A|2^'^'^~^\ consistently with the bounds 
on derived in the previous section. 

The bound Eg. (14. 5 p is one of the basic ingredients from which the main 
results of this paper follow. It makes apparent that all the terms with at 
least one endpoint on scale or with one endpoint of type u are strongly 
suppressed, in the sense that they are proportional to a short memory factor 
2e(h-N)^ which go to zero in the limit — )■ oo (which corresponds to the 
scaling limit a — 0): these are the terms leading to the correction term 
i?(")(xi,ji; ■ ■ ■ ;Xm,jm) in Theorem 11.2^ as explained in more detail below. 
The dominant terms are those coming from trees r G T^'^om' which are 
of order 1 with respect to A (because the size of the special endpoints is 
measured by which is analytically close to 1), and whose dimensional 
bound does not vanish in the N ^ oo limit. In the following, we show that: 
(1) the contributions to the m point function from trees r G TJ^^o m can be 
written as the expression in the second line of Eg. (II. 15]) plus a remainder, 
bounded by the r.h.s. of Eg. (II. 16]) : (2) the sum of the contributions from all 
the trees with n > 1 or uiq < m can be bounded by the r.h.s. of Eg. (II. 161) . 
This will conclude the proof of the main results stated in the introduction. 

4.1 The dominant contributions 

As discussed above, the dominant contribution to the m point energy cor- 
relation function come from trees r G T^'^dm- Among these, we further 

distinguish the contributions from trees in T^^l ^ with at least one endpoint 
on scale N + 2, which will be discussed at the end of this subsection, and 
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those with all the endpoints on scales < + 2, which can be written as 



N 



if(xi, . . . ,x^) := m! ^ ^ S'^,p(xi, ji; ■ ■ ■ ;x^, j^) , 

h=h„-l ^pr^h) FgVt: 

(4.6) 

where the * on the sum indicates the constraint that all the endpoints are 
on scale < + 2. By summing over trees, and by explicitly computing the 
tree values, K{'Ki, . . . , x^) can be written in a way very similar to Eg. (11.81) : 

1 / i \ 

K(xi,...,x„) = -— - • (4.7) 



2m Vtt 



'ren^i,...,^}'^!.-.^^ j=i 



where the sum over the h/s runs over < hj < N, with the convention 
that g^'^'^'^ is the propagator of Equation (14. 7p differs from the corre- 

sponding expression for the free system in the renormalization of the special 
endpoints, -^^j^nf/i ft +i}' ^^'^ propagators. 

First of all, using Eq. (l3.6ip . we can rewrite h +i} ~ ^-LW + 

0(|A|2^^). Similarly, using the explicit expression for c/W, Eq.([336D, and the 
bounds Eqs. (l3.55p . (l3.6ip . we rewrite (after having taken the thermodynamic 
limit L — > oo) g^'^\x — y) = g'^''^(x — y) + r/i(x — y) where, for all ha- < h < N, 

^ ^' ^-oo(A) 27r|z)^(k)|2+[z*(A)a]H-^^*(^V ^a(k) 

(4.8) 

where rh is a correction bounded as 



x - y)|| < C|A|2^(^-^)2V'=2'l"-yl . (4.9) 



li h = h„ we still write g^^'^\'x) = q^^'^\x) + r/i^(x), with r^^ bounded as in 
Eq.dM]) and g(''-)(x) defined by an expression analogous to Eq. (l4.8p . with 
the only difference that //i^(k) is replaced by J2k<h fk(X) = ■ 

Correspondingly, we rewrite Eq. (l4.7p as the term obtained by replacing 

all the factors Z^in{hi h^+i} ^-loW and all the propagators g^^^ by g''^^ 
plus a remainder, which is the sum of all the terms involving at least one 
factor Zj^^ — Z^l^, or one factor r/i(x): 

K{xi, . . . , x^) = K^'\xu ...,xm) + K'ixi, . . . , x^) (4.10) 
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(a) 




h K hi h2 hs /i4 N + 2 

(b) 

Figure 9: (a) A diagram representing one of the terms in the expression in 
Equation (14. 6p for il'(xi, . . . , xg), with the scale labels of the propagators 
shown. The vertices are grouped into clusters according to scales of the 
connecting propagators, (b) A GN tree r associated with this term; note 
that the branching points of the trees correspond to the clusters of vertices 
shown above, and that this tree is uniquely specified apart from the choice 
of which endpoint is identified with which coordinate. 
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where 



irW(x„...,xJ = -^(l)". (4.11) 



m 



'rGn{i,.,.,„} t-"!.-.^™ j = l 

ni,...,/lm 

Using the fact that J2h<N fhi.^) = 1 we see that Eg. (14. Ill) is equal to the 
expression in the second hne of Eq. Ol.lSp . provided that Z{X) in Eq. (ll.l4|) is 
fixed to be Z{X) = zi^i(A)/Z_oo(A). Moreover, using the bounds Eqs. (l3.37p . 
fl3.6ip . and (14.91) on g'^'^\ zj^^ — Z^l^, and r^, respectively, the correction term 
K'{yii, . . . , Xm) can be bounded as: 

|ir'(xi,...,x^)| < C""|A|m! J] 2^('^"'"-^)[ JJ 2(^^-^)h.e~''''^^^^ , 

(4.12) 

where the sum over trees has been re-introduced in order to properly keep 
track of the hierarchical scale structure induced by the choice of the scale 
labels of the propagators, as explained in the caption to Fig|9l Moreover 
Vm is the rightmost non trivial vertex of r (e.g., = in the example of 
FigJH]) and the set Vnt{T) is the set of non trivial vertices of r. Note that 
^v€V„t{T) '^^"^v > 2''o£), where := h^* and D = is the diameter of 
X = {xi, . . . , Xm}. Therefore, 

|ir'(xi,...,x^)| < C™|A|m!^2" ^ ^-c'2^oD^e(h.^-N) _ 



' N:0,m 



[ 



Yl 2(^-i)'^''e-^'2'"^''J , (4.13) 

where c' can be chosen to be half of the constant c in Eq. (l4.12l) . We then 
separate the sum over r into a sum over the various scale labels, and a sum 
over the remaining structure of the tree, which is indexed by an unlabeled 
tree t having the same endpoints as r. We denote by To,m the set of un- 
labeled trees with n = normal endpoints and m special endpoints. Then 
introducing the variable j = h — Hq, ignoring all the constraints in the sum 
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over scales other than h < h^, 



oo oo 
hvm=~'°° veVntir): hv=—oo 

Now, the sum over j gives a constant and we perform the other sums as 
follows. Letting —hs^ be the integer part of log2 so that 2~^^^ < Sy < 
2(i-'^«J, then for any > 0, 



hv=—oo hv=~oo 

(4.14) 

where 



y 2"V^'2'^< rdx2"^e-^'2>maxy"e-'='^ = (c')-4^^+f-) 



/i=— c» 



(4.15) 

In particTilar, if o; > 1, then < C"q;" for a suitable C > 0. Therefore, 
letting (5 = 5x = mini<j<j<TO |xj — Xj\ and using also the fact that > 

{Sv - 1)(^, 

|X'(xi,...,x„)|<(C")'"|A|m!2-^^ J] ( 



Recall that s^* > 2: under this condition it is easy to check that 



Finally, recalling that X]i.ey„t(t)('^f ~ ^) = rn — 1 and noting that the sum 
over t simply gives the number of unlabeled trees with m endpoints (which 
is no more than 4™), we get 

|X'(xi, . . . , xJI < (C")"^|A|m!^ (^)'(-^)' , (4.17) 
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which is the desired estimate on the contributions to the m points function 
coming from the trees in Tj^^g ^. 

The contributions to the m point function coming from trees r G Tj^^o m 
with at least one special endpoint on scale N + 2 can be bounded by a similar 
expression, with the difference that the factor |A| in the r.h.s. of Eq. fl4.17p 
should be replaced by 1. We do not belabor the details of this computation, 
which is completely analogous to the one leading to Eq.f l4.17p . Let us just 
remark that these are the terms responsible of the fact that the r.h.s. of 
Eq. fll.l6p is of order 1 w.r.t. A. Of course, if desired, these contributions 
can be written explicitly (possibly modulo further corrections bounded as 
in Eq.f l4.17]) ). by making use of the explicit expression of the source term 
B{ip, X, A), which is obtained from Eq. fl2.45p by performing the linear change 
of variables Eq. f l2.52p , and of the explicit expression of the propagator of the 
X field, Eq. (12.48]) . If we decided to isolate these terms from the correction 
then the remaining contributions would be bounded as in Eq. (l4.17]) . 
as it follows from Eq. (l4.17p itself and from the discussion in the following 
subsection. 



4.2 The subdominant contributions 

It remains to bound the trees with n > 1, 

K"{xi,ji;--- ;x^,jrn) ■■= (4.18) 

N 

■=m\ ^ ^ 5Z '5r,p(xi,ji;--- ;x^,j^) , 

h=ha-l n>l ^pT-C") PG-Pr: 

mo<m N:n,mo |p^^|^|pA|^^ 

where S't-_p(xi, ji; ■ ■ ■ ; x^, jm) can be estimated as in Eq. (l4.5p . For simplicity, 
we start by looking at the contributions to K" coming from trees with m = 
mo, to be called K'": 

i^:'"(xi,ji;-- - ;xm,j™) := (4.19) 

N 

■=m\ Y Y Yl Yl '^r,p(xi,ii;--- ;x„,j„) , 
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which is bounded as: 

N 

\K"'{^^,j,;---;^m,Jm)\<m\ E 

h=h^-l n>l ^^(h) PGPr: 
_2hi2-ni)20{h+-N)^ JJ' (-(^1 | ■)max{l,c|P„* | } _ (4.20) 

feV„t(r*) vnote.p. 

where is the scale of the rightmost normal endpoint and En{T) is the 
set of normal endpoints of r. We now want to reduce this expression to a 
simplified form as close as possible to Eq. fl4.12p and then bound the simplified 
expression by a strategy similar to the one used in the previous subsection. 
The first step consists in "pruning" the tree r of the branches that are not 
in r*. To this purpose, we make the following rearrangement. Let T^h*-m 
the set of labelled trees with m special endpoints, with no normal endpoints, 
and whose leftmost vertex Vq is nontrivial and on scale h^. If r* G h*-m 

and h < /ig, let Tl^\T*;n) C T^'^lj^ be the set of trees r with root scale h 
and n normal endpoints such that r* is the subtree of r connecting its special 
endpoints. Fix any 6 e (0, 1), e e (0, 1/2) and (p e (max{0, 1 - 6^ - Se}, 1 - 6') 
and define, recalling that m„ = \P^\, 

{1 + £ — m^, > 1 
-e-(t), = 1 . (4.21) 

0, m.u = 

It is easy to check that, if \Pf \ > 2, and recalling that d^i^Pv) = min{l — 
\\Pf\-m,,-l}, then 

- I = d,{P,) - UP:^) < -1(1 - e - <f))\PS\ - 95^^,0 . (4.22) 

o 

It is convenient to rewrite Eq. fl4.20p as 

N 

|ir'"(xi,ji;---;x„,j^)|<C'"m! 2^^^'-^^ ^ ■ 

h=h„-l h<h*<N+l T*eT* 

. ^ 2eih+~N) ^ J j-j- 2K-rfO('i--V)j J Y[ (C|A|)'^^^-ti''=l-P-'^l> 

-Tug | — KUQ |— ''i 

66 




where V{t) is the set of vertices of r that are not endpoints, and similarly 
for V{t*). Using Eq. fl4.22l) we can perform the sum over r G T^^\t*; n) and 
over P G Pr, which gives 

where is the highest scale in r* . Plugging Eq. fl4.23p back into the previous 
bound on K'" , and summing over n > 1 and over h < gives: 

■[ n 2^^''*"-^^''''e"''^'"'''''j 1^ JJ^ 2'^'''(''"~V) _ (4.24) 

feV„t(r*) t)ey(r*) 

We now go from the sum over r* to an estimate in terms of a "contracted 
tree", called r**, which is obtained from r* by removing the set E of vertices 
which precede exactly one special endpoint, see Figure [TUl The sum over r* 
can then be expressed as the sum over r** (trees with only special endpoints, 
where each endpoint is attached to a branching point) and the sum over 
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compatible E (i.e. over the number of vertices inserted before each special 
endpoint); then Eq. fl4.24p becomes (calling T^*/j*.„ the set of contracted trees 
r** with root on scale hQ — 1 and m endpoints) 



l^'"(xi,ji; 



) -^m 1 Jm) 



<C"^|A|m! 2^*' 



(2-m) 



E 



T"er* 



fGV„t(T") 



dv(hv—h^i) 



vGV{t") 

2e{hm-N) 



(4.25) 



v€E 



We use the factors 2'^^^^^'^^'^ to perform the summation over E, and use the 
remaining factor to replace hm with the scale of the last branching point of 
T**, which we denote by h^, see FiglTOl 



5Z 1^ Jj2~(''+'^)(''""^'''^ 

This leaves us with the estimate: 

|ir'"(xi,ji;---;x™,j„)|<C"^|A|m! J] 2'^S(2-H 

?t*<Ar+i T"er; 

2^'(^"-JV) 



(4.26) 



E 



1^ 2(iv{hv-h^,) 

vev{T**) 



(4.27) 



which is similar to Eq. fl4.12p . We would be tempted to proceed as we did 
after Eq. fl4.12p . that is by ignoring all the constraints in the sum over the 
scale labels and reduce the r.h.s. of Eq. (14.271) to a product of terms of the 
form Yl'h=-oo 2""''*"e~'^^ , which is fine provided that > 0. However, here 
the exponent ay associated with Vq is equal to 2s** — m, which in general is 
negative. 

To rectify this we engage in a further rearrangement of the estimate, 
for which we introduce the following notation. Let S*'"^ denote, for v G 
Vnt{T**), the set of branching points following v but not following any vertex 
w G Vnt{T**) such that w > v (see Figure [TT]), S*'^ the set of endpoints 
immediately following v, and S* = S*'^US*''^; also let s* = \S*\ and similarly 
for s*'^ and s*'^. Since we have = l + e — m^ for all non-endpoint vertices, 
we have the following estimate for each v G S*'?', 



(4.28) 
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Figure 11: Example of the sets S*'^ and S*'"^. 



which since y^,,,^Q*,2(mv — l — e) — m — s*'» —(l + e)s*l = m—(l + e)s*,*+es*,l 
imphes 

^g5*;2 V*<W<V 



"0 



"0 

We now have a factor of 2(^*''~"*''~-'^+^)'*'' for each v G S*l^, and this can still 
produce a divergence. For these vertices we have 



which we apply inductively to all branching vertices to obtain 



|/r'"(xi,ji;...;x^,jJ|<C"^|A|m!2-^^ ^ [ J] 2 



/i5<Af+l DeV„t(r") 

(4.31) 



* * ^ '7"'* * 
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with 



till 



{l-e)sl + esl'\ v = v* 

(1 - e){sl - 1) + £s*/, otherwise 



(4.32) 



Note that we then have > 1/2 and ^ = m + ^. We can now sum this 
in the same way as explained after Eg. (14. 121) : noting that ay* > 2 — 2e, this 
gives 



which is vahd for any 6 G (0, 1) and any e G (0, 1/2), but that the "constants" 
depend on 9, e, so we cannot take 6* — )■ 1 or e — )■ 0. 

The contributions to the m points energy correlations from trees with 
n > 1 normal endpoints and mo < m special endpoints can be estimated via 
an analogous procedure and we are lead again to a bound like Eq.f l4.33p . We 
do not belabor the details of this computation here. Therefore, the proofs of 
Theorems 11.11 and 11.21 are complete. 
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